ESD-TR-87-028 

ESL-TR-85-039 


Final  Report 


Laser  Remote  Sensing 
of  Atmospheric  Pollutants 


30  September  1985 


Lincoln  Laboratory 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 

Lexington,  Massachusetts 


Prepared  for  the  Air  Force 

under  Electronic  Systems  Division  Contract  F19628-85-C-0002. 


Approved  for  public  release;  distribution  unlimited. 


The  work  reported  in  this  document  was  performed  at  Lincoln  Laboratory,  a  center 
for  research  operated  by  Massachusetts  Institute  of  Technology.  This  work  was 
supported  by  the  Department  of  the  Air  Force  under  Contract  F19628-85-C-0002,  in 
part  with  specific  funding  from  the  Engineering  and  Services  Laboratory,  Air  Force 
Engineering  and  Services  Center,  Tyndall  Air  Force  Base,  Florida. 

This  report  may  be  reproduced  to  satisfy  needs  of  U.S.  Government  agencies. 


The  views  and  conclusions  contained  in  this  document  are  those  of  the  contractor 
and  should  not  be  interpreted  as  necessarily  representing  the  official  policies, 
either  expressed  or  implied,  of  the  United  States  Government. 


The  ESD  Public  Affairs  Office  has  reviewed  this  report,  and 
it  is  releasable  to  the  National  Technical  Information 
Service,  where  it  will  be  available  to  the  general  public, 
including  foreign  nationals. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 
FOR  THE  COMMANDER 


Thomas  J.  Alpert,  Major,  USAF 

Chief,  ESD  Lincoln  Laboratory  Project  Office 


Non-Lincoln  Recipients 

PLEASE  DO  NOT  RETURN 

Permission  is  given  to  destroy  this  document 
when  it  is  no  longer  needed. 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
LINCOLN  LABORATORY 


LASER  REMOTE  SENSING 
OF  ATMOSPHERIC  POLLUTANTS 


D.K.  KILLING ER 
N.  MENYUK 
A.  MOORADIAN 
Group  82 


FINAL  REPORT 
TO  THE 

AIR  FORCE  ENGINEERING  AND  SERVICES  CENTER 


1  OCTOBER  1984  —  30  SEPTEMBER  1985 

ISSUED  21  JULY  1987 


Approved  for  public  release;  distribution  unlimited. 


LEXINGTON 


MASSACHUSETTS 


ABSTRACT 


This  FY85  Final  Report  on  the  program  entitled  “Laser  Remote 
Sensing  of  Atmospheric  Pollutants”  covers  the  period  1  October  1984 
to  30  September  1985.  The  work  was  supported  by  the  Air  Force 
Engineering  and  Services  Center. 


TABLE  OF  CONTENTS 


Abstract  iii 

List  of  Illustrations  vii 

I.  INTRODUCTION  1 

II.  Co:MgF2  DIAL  Measurements  of  Atmospheric  Species  1 

III.  Remote  Sensing  of  HC1  Plume  4 

IV.  Comparison  of  Experimental  Data  with  Predictions  5 

V.  Associated  Programs  in  Co:MgF2  Laser  Development  6 

VI.  Conclusion  7 

Appendix  A:  Paper  entitled  “Atmospheric  Remote  Sensing  of  Water  Vapor, 

HC1  and  CH4  Using  a  Continuously  Tunable  Co:MgF2  Laser”  9 

Appendix  B:  Paper  entitled  “Frequency  Tuning  Characteristics  of  a 

Q-switched  Co:MgF2Laser”  27 

Appendix  C:  Paper  entitled  “Laser  Remote  Sensing  of  the  Atmosphere”  37 


v 


LIST  OF  ILLUSTRATIONS 


Figure 

No. 

1 

2 

3 

4 

5 

6 


Page 


Photograph  of  Co:MgF2  Laser  2 

Photograph  of  Co:MgF2  DIAL  System  3 

Range-Resolved  Co:MgF2  Lidar  Returns  Backscattered  from 

Atmospheric  Aerosols  as  a  Function  of  Range  4 

Differential-Absorption  (Range-Resolved)  Remote  Sensing  of  HC1 

Plume.  The  Lidar  Beam  Passed  Through  a  Semi-Enclosed  Cell  at  a 

Range  of  200  m  Which  Contained  the  Source  of  the  HC1  Plume  4 

Predicted  Detection  Range  for  Path- Averaged  Remote  Sensing  of  HC1 

in  the  Atmosphere  5 

Predicted  Detection  Range  for  Range-Resolved  Remote  Sensing  of 

HC1  in  the  Atmosphere  6 


Vll 


LASER  REMOTE  SENSING 
OF  ATMOSPHERIC  POLLUTANTS 


I.  INTRODUCTION 

This  is  the  FY85  Final  Report  on  the  program  entitled  “Laser  Remote  Sensing  of 
Atmospheric  Pollutants”  supported  by  the  Air  Force  Engineering  and  Services  Center. 

The  specific  tasks  which  were  conducted  during  FY85  for  this  research  program  consisted 
of  the  following: 

(a)  The  continuation  of  the  Co:MgF2  DIAL  measurement  of  atmospheric  species, 
and 

(b)  Initial  remote  sensing  demonstration  of  a  HC1  plume  in  the  atmosphere. 

Each  of  these  tasks  is  described  in  detail  in  the  following  sections.  Supportive  documentation 
is  included  in  the  appendices. 

II.  Co:MgF2  DIAL  MEASUREMENTS  OF  ATMOSPHERIC  SPECIES 

The  Co:MgF2  laser  differential-absorption  lidar  (DIAL)  system  reported  in  the  FY84  Final 
Report*  was  modified  in  order  to  increase  laser  output  power  and  to  incorporate  a  better  tele¬ 
scope;  Figures  1  and  2  show  photographs  of  the  Co:MgF2  laser  and  DIAL  system.  These  modifi¬ 
cations  enhanced  the  lidar  signal  sufficiently  so  that  returns  from  atmospheric  aerosols  were  eas¬ 
ily  observable.  As  a  result,  range-resolved  laser  remote  sensing  measurements  of  the  atmosphere 
were  made  using  this  system  and  were  used  to  establish  detection  range  and  sensitivity.  The  tech¬ 
nical  details  of  the  laser  system  have  been  presented  in  the  FY84  Final  Report,  and  details  of  the 
DIAL  experiments  and  related  investigations  are  given  in  Appendices  A  and  B.  A  synopsis  of  the 
study  is  given  below. 

The  Co:MgF2  DIAL  system  had  sufficient  laser  power  (10  mJ/pulse)  and  optical  detection 
sensitivity  to  permit  the  remote  sensing  of  lidar  backscatter  from  aerosols  and  molecules  in  the 
atmosphere.  Figure  3  shows  a  plot  of  the  Co:MgF2  lidar  returns  from  atmospheric  aerosols  as  a 
function  of  range.  As  seen  in  the  figure,  lidar  returns  were  recorded  at  ranges  out  to  3.8  km, 
with  significant  signal-to-noise  ratios  measured  at  ranges  out  to  about  3  km.  These  results  are  in 
good  agreement  with  theoretical  calculations  which  predicted  that  the  Co:MgF2  DIAL  system 
should  be  capable  of  detecting  atmospheric  backscatter  returns  out  to  ranges  of  ~3  to  4  km  (see 
Appendix  A). 


*  D.K.  Killinger,  N.  Menyuk,  and  A.  Mooradian,  “Laser  Remote  Sensing  of  the  Atmosphere,” 
FY84  Final  Report,  Lincoln  Laboratory,  MIT,  ESD-TR-85-212,  ESL-TR-85-34  (30  September 
1984),  DTIC  AD-A159531/3. 
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Figure  1.  Photograph  of  Co:MgF2  laser. 
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Figure  2.  Photograph  of  Co:MgF2  DIAL  system. 
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Figure  3.  Range-resolved  Co:MgF2  lidar  returns  backscattered  from  atmospheric  aerosols 
as  a  function  of  range- 


ill.  REMOTE  SENSING  OF  HCI  PLUME 

The  Co:MgF2  DIAL  system,  operating  near  1.7525  jum,  was  used  for  the  range-resolved 
remote  sensing  of  HCI  in  the  atmosphere.  Previous  DIAL  experiments  (FY84  Final  Report)  had 
established  that  a  negligible  amount  of  ambient  HCI  (<0.2  ppm)  was  present  in  the  atmosphere. 
Therefore,  an  artificial  source  of  HCI  was  used  in  order  to  provide  a  plume  of  HCI  suitable  for 
remote  sensing.  In  order  to  reduce  the  amount  of  gaseous  HCI  injected  into  the  atmosphere,  a 
semi-enclosed  cell  was  used.  The  1-m-long,  0.8-m-diam.,  windowless  cell  was  placed  at  a  distance 
of  200  m  from  the  DIAL  system,  and  the  0.3-m-diam.  lidar  beam  passed  unobstructed  through 
the  cell.  Figure  4  shows  the  normalized  lidar  signal  obtained  after  HCI  was  injected  into  the  cell. 


0  200  400  600  800  1000 
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Figure  4.  Differential-absorption  (range-resolved)  remote  sensing  of  HCI  plume.  The  lidar  beam  passed 
through  a  semi-enclosed  cell  at  a  range  of  200  m  which  contained  the  source  of  the  HCI  plume. 
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The  difference  between  the  on-line  and  off-line  signals  indicates  that  ~40  ppm  of  HC1  was  pres¬ 
ent.  The  minimum  differential-absorption  signal  observable  was  a  few  percent  of  this  value.  As 
such,  we  estimate  that  the  sensitivity  of  the  Co:MgF2  DIAL  system  was  ~1  to  3  ppm  with  a 
200-m  range  resolution.  The  range  resolution  is  determined  by  the  laser  pulse  length  and,  as  seen 
in  the  figure,  is  responsible  for  the  considerable  degree  of  spatial  averaging  seen  in  the  data.  A 
shorter  pulse  length  of  the  Co:MgF2  laser  would  permit  a  shorter  range  resolution  tp  be 
obtained. 

IV.  COMPARISON  OF  EXPERIMENTAL  DATA  WITH  PREDICTIONS 

One  can  compare  the  experimental  results  shown  in  Figure  4  with  the  theoretical  predictions 
of  the  sensitivity  of  the  Co:MgF2  DIAL  system.  The  predictions  are  presented  in  Figures  5  and  6 
which  show  the  calculated  minimum  detectable  concentration  of  atmospheric  HC1  as  a  function 
of  range  for  both  path-averaged  and  range-resolved  measurements  assuming  a  laser  pulse  energy 
of  10  mJ;  details  of  the  calculations  are  given  in  Appendix  A. 

Figure  5  indicates  that  the  system  should  be  capable  of  detecting  path-averaged  concentra¬ 
tions  of  HC1  between  0.1  and  1  ppm  over  ranges  to  about  10  km.  This  prediction  is  in  good 


Figure  5 .  Predicted  detection  range  for  path-averaged  remote  sensing  of  HCl  in  the  atmosphere. 


5 


Figure  6.  Predicted  detection  range  for  range-resolved  remote  sensing  of  HCl  in  the  atmosphere. 


agreement  with  the  experimental  data  which  indicated  a  sensitivity  on  the  order  of  0.2  ppm  and 
detection  ranges  out  to  6  km  (see  Appendix  A). 

The  range-resolved  aerosol  backscatter  predictions  shown  in  Figure  6  indicate  that  a  concen¬ 
tration  sensitivity  of  a  few  pulses  per  million  to  a  range  of  ~2  km  may  reasonably  be  expected. 
These  predictions  are  in  good  agreement  with  the  experimental  data  shown  in  Figure  4. 

V.  ASSOCIATED  PROGRAMS  IN  Co:MgF2  LASER  DEVELOPMENT 

The  detection  range  of  the  Co:MgF2  DIAL  system,  while  certainly  adequate  for  many 
remote  sensing  applications,  could  be  enhanced  by  increasing  the  output  power  of  the  Co:MgF2 
laser.  The  current  laser  is  limited  in  energy  to  ~10  mJ/ pulse  because  of  optical  damage  to  the 
laser  cavity  mirrors;  at  higher  Q-switched  powers,  the  dichroic  mirror  is  irreversibly  damaged. 
Several  studies  have  been  made  by  others  to  increase  the  power  out  of  the  Co:MgF2  laser. 
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Sanders  Associates*  has  recently  completed  a  program  funded  by  AFSC/ Wright  Patterson 
Avionics  Laboratory  to  develop  a  high-power  Co:MgF2  laser.  This  program  recently  demon¬ 
strated  30-Hz  pulsed  laser  operation  at  the  peak  of  the  Co:MgF2  gain  (1.9  jum)  of  0.5  J/pulse  in 
the  long-pulsed  (non-Q-switched)  mode  and  150  m J/pulse  in  a  150-ns  Q-switched  pulse.  Higher 
powers  were  limited  by  optical  damage  to  the  cavity  mirrors. 

Finally,  it  should  be  noted  that  the  Co:MgF2  laser  is  currently  being  investigated  also  by 
Schwartz  Electro-Opticst  with  funding  from  NASA/ Langley.  This  program  has  demonstrated 
operations  of  the  Co:MgF2  laser  at  0°C  using  a  thermoelectric  cooler  instead  of  liquid  nitrogen 
cooling  for  the  Co:MgF2  crystal. 

VI.  CONCLUSION 

The  Co:MgF2  DIAL  system  has  been  shown  to  be  useful  for  path-averaged  and  range- 
resolved  remote  sensing  of  atmospheric  species  at  ranges  up  to  6  and  3  km,  respectively.  Range- 
resolved  DIAL  measurements  of  HC1  plumes  have  established  a  detection  sensitivity  of  a  few 
parts  per  million  and  indicate  that  the  Co:MgF2  DIAL  system  can  be  useful  for  the  detection  of 
a  wide  variety  of  HC1  gaseous  sources  in  the  atmosphere.  Future  developments  related  to  enhanc¬ 
ing  the  output  power  of  the  Co:MgF2  laser  should  certainly  increase  the  detection  range  and  sen¬ 
sitivity  of  the  DIAL  system. 


*  M.  Thomas,  Sanders  Associates  (private  communication). 

+  P.  Moulton,  Schwartz  Electro-Optics  (private  communication). 
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APPENDIX  A 


The  following  is  a  reprint  of  a  journal  article  accepted  for  publication  in  Applied  Optics, 
entitled  “Atmospheric  Remote  Sensing  of  Water  Vapor,  HC1  and  CH4  Using  a  Continuously 
Tunable  Co:MgF2  Laser.” 


ATMOSPHERIC  REMOTE  SENSING  OF  WATER  VAPOR,  HC1  AND  CH4 
USING  A  CONTINUOUSLY  TUNABLE  Co:MgF2  LASER 


N.  Menyuk  and  D.K.  Killinger 
Lincoln  Laboratory,  Massachusetts  Institute  of  Technology 

Lexington,  MA  02173-0073 


ABSTRACT 


A  differential-absorption  lidar  (DIAL)  system  has  been  developed  which 
uses  a  continuously  tunable  (1.5  —  2.3  jum)  cobalt-doped  magnesium  fluoride 
(Co:MgF2)  laser  as  the  radiation  source.  Preliminary  atmospheric  measure¬ 
ments  of  water  vapor,  HC1  and  CH4  have  been  made  with  this  system, 
including  both  path-averaged  and  range-resolved  DIAL  measurements  at 
ranges  up  to  6  km  and  3  km,  respectively. 


*  This  work  was  supported  by  the  National  Aeronautics  and  Space  Administration 
and  the  Air  Force  Engineering  and  Services  Center. 


H 
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I.  INTRODUCTION 


There  is  an  acknowledged  need  for  tunable  high-power  lasers  in  the  near-to  middle-IR  to 
serve  as  sources  in  lidar  systems  for  remote  sensing  of  the  atmosphere.1  Recently  several  new  IR 
lasers  which  show  potential  for  use  as  sources  for  lidar  or  differential-absorption  lidar  (DIAL) 
systems  have  been  developed.  They  include  optical  parameteric  oscillators  (OPO)  pumped  by  a 
Nd:YAG  laser,2  Ti:Al203  lasers,3  and  OH:NaCl  lasers.4 

The  solid-state,  cobalt-doped  magnesium  fluoride  laser  (Coih4gF2),5>^  which  is  continuously 
tunable  from  1.5  to  2.3  /im,  also  falls  in  this  category.  In  this  paper  we  report  on  the  adaptation 
and  use  of  a  Co:MgF2  laser  as  a  practical  pulsed  source  of  coherent  radiation  in  a  DIAL  system 
and  describe  the  use  of  this  system  to  measure  atmospheric  species  in  the  atmosphere.  The  system 
can  produce  either  path-averaged  or  range-resolved  remote  sensing  measurements  to  ranges  of 
6  km  and  3  km,  respectively. 

II.  Co:MgF2  DIAL  SYSTEM 

The  DIAL  system  is  described  in  this  section.  In  addition  to  the  Co:MgF2  laser,  the  system 
includes  optical  diagnostics  for  energy  and  spectral  monitoring  of  the  laser,  lidar  optics  consisting 
of  a  beam  steering  mirror,  receiving  telescope,  and  optical  detector,  and  a  computer  to  record  the 
DIAL  signals  and  to  control  and  monitor  the  laser  wavelength.  A  schematic  of  the  complete  sys¬ 
tem  is  given  in  Figure  1. 

A.  Co:MgF2  Laser 

Early  investigations  of  the  Co:MgF2  laser  by  Moulton5  established  that  the  laser  was  con¬ 
tinuously  tunable  from  1.5  to  2.3  /im,  could  operate  either  cw  or  pulsed,  and  could  radiate  at  a 
single  frequency.  Lovold  et  al.6  modified  the  original  design  so  that  the  laser  could  be  Q-switched 
and  accurately  controlled  in  wavelength  in  order  to  make  the  Co;MgF2  laser  a  practical  radiation 
source  for  remote  sensing  measurements.  The  resulting  configuration  is  shown  in  the  upper  left 
portion  of  Figure  1. 

The  Co:MgF2  laser  was  longitudinally  pumped  by  a  1.32  /im  Nd:YAG  laser  which  emitted 
400  mJ,  800  /xs  duration,  linearly  polarized,  multimode  pulses  at  a  3  Hz  pulse  repetition  fre¬ 
quency.  The  pump  pulse  was  focused  onto  a  27  mm  long  Co:MgF2  crystal  cut  at  the  Brewster 
angle.  The  crystal  was  mounted  in  a  liquid  nitrogen  Dewar  which,  in  turn,  had  entrance  and  exit 
windows  set  at  the  Brewster  angle.  The  output  coupling  mirror  (M2  in  Figure  1)  used  for  the 
remote  sensing  measurements  reported  here  had  a  95%  reflectivity  at  1.7  Mm,  and  the  dichroic 
back  mirror  (Ml  in  Figure  2)  had  high  reflectivity  at  1.7  Mm  and  high  transmission  at  1.32  Mm. 
An  intracavity,  fused-silica,  acousto-optic  Q-switch  was  used  to  reduce  the  laser  pulse  length  to 
300  ns;  when  operated  in  a  free-running,  pulsed  mode  (non-Q-switched),  the  output  appears  as  a 
series  of  spikes  spaced  approximately  50  ms  apart  and  covering  a  time  span  of  about  300  ms, 
making  that  operative  mode  difficult  to  use  for  remote  sensing  application. 
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Course  tuning  and  linewidth  control  of  the  the  laser  was  obtained  using  a  three-element 
birefringent  tuning  filter6  which  reduced  the  linewidth  to  approximately  1  cm'1  and  permitted 
laser  tuning  over  the  spectral  range  from  1.6  to  1.9  /i m,  the  maximum  range  obtainable  with  a 
single  set  of  laser  cavity  mirrors.  A  quartz  intracavity  etalon,  0.25  mm  thick  and  coated  for  20% 
reflectivity,  reduced  the  linewidth  to  0.15  cm-1  and  was  angle  tuned  for  precision  scanning  of  the 
laser  wavelength  over  a  limited  range  of  the  order  of  5  cm'1.  The  output  energy  of  the  Co:MgF2 
laser  was  limited  to  approximately  10  mJ  when  operated  Q-switched  in  order  to  avoid  optical 
damage  to  the  dichroic  mirror. 

B.  DIAL  System 

As  seen  in  Figure  1,  the  bulk  of  the  Co:MgF2  laser  output  is  used  for  remote  sensing,  but 
portions  of  the  output  are  sampled  using  beamsplitters  to  provide  diagnostics.  In  particular,  a 
room  temperature  HgCdTe  detector7  monitors  the  energy  of  each  pulse  and  provides  an  optical 
trigger  to  various  components,  and  the  1-m  long  absorption  cell  shown  in  the  figure  permits 
laboratory  calibration  of  the  absorption  of  the  molecule  being  investigated  as  a  function  of  the 
laser  wavelength.  The  spectrometer  shown  monitors  the  laser  frequency  at  a  particular  position; 
however,  when  the  laser  is  used  in  a  continous  frequency  scanning  mode,  the  fixed  air-gap  Fabry- 
Perot  gives  an  accurate  measure  of  the  laser  wavelength  variation,  which  is  fed  directly  to  the 
computer.  In  addition,  the  spectrometer  is  used  as  a  spectrograph,  with  the  output  imaged  on  a 
64-element  pyroelectric  detector  array  with  100  /xm  spacing  of  the  individual  elements.  This 
arrangement  permits  visual  observation  in  real  time  of  the  laser  wavelength  output  over  a  range 
of  about  25  cm-1  to  ensure  that  the  output  is  a  single  mode  of  the  intracavity  etalon  and  is  con¬ 
tinuously  tuning  over  the  free  spectral  range  of  the  etalon. 

The  remote  sensing  portion  of  the  output  beam  is  expanded  to  a  3  cm  diameter  and  is 
directed  to  the  target  with  a  50  cm  beam  steering  mirror.  The  return  lidar  beam  is  sent  through  a 
30  cm  Cassegrain  telescope  and  detected  by  a  cooled  InSb  detector.  The  output  signal  is  then 
measured  by  an  A/D  converter  or  a  waveform  digitizer  and  analyzed  by  the  computer.  The  com¬ 
puter  also  controls  and  scans  the  wavelength  of  the  laser  by  use  of  a  Burleigh  “Inchworm”  which 
changes  the  slant  angle  of  the  thin  etalon  in  the  laser  cavity. 

The  system  can  run  the  wavelength  control  of  the  laser  either  in  a  continuous  scanning 
mode,  or  in  a  differential-absorption  (DIAL)  mode.  In  this  latter  case  only  two  wavelengths,  on 
and  off  an  absorption  line,  are  employed  and  the  position  of  the  etalon  is  toggled  between  these 
two  positions.  In  the  continuous  scanning  mode,  the  output  of  the  Fabry-Perot  is  used  to  cali¬ 
brate  the  relative  wavelength  of  the  laser  and  to  correctly  scale  the  wavelength  scan  of  the  laser 
and  lidar  returns. 

III.  PATH-AVERAGED  DIAL  MEASUREMENTS 

The  system  has  been  operated  in  the  continuous  scanning  mode  to  obtain  path-averaged 
(column  content)  measurements  using  backscattered  returns  from  topographic  targets.  Figure  2 
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shows  a  plot  of  the  normalized  lidar  returns  obtained  from  foliage  on  hillsides  located  3  km  and 
6.7  km  from  the  laboratory  as  a  function  of  the  laser  wavelength.  The  wavelength  scan  covers 
approximately  5  cm*1.  All  the  major  absorption  regions  can  be  identified  with  known  water 
vapor  absorption  lines.8"10  While  more  than  one  absorption  line  is  involved  in  most  of  the  transi¬ 
tion  regions,  the  absorption  near  1.7678  pm  is  primarily  due  to  a  single  H20  transition  with  an 
absorption  linewidth  equal  to  0.11  cm"1  at  atmospheric  pressure  and  an  expected  peak  absorption 
coefficient,  aa,  of  5.2  X  104  (atm-cm)"1.8  One  deduces  from  the  observed  absorption  strength  of 
52%  in  Figure  2(a)  that  approximately  0.0023  atmospheres  (1.8  Torr)  of  water  vapor  was  present 
in  the  atmosphere. 

This  DIAL  measurement  can  be  compared  to  that  obtained  using  a  wet-bulb/ dry-bulb  ther¬ 
mometer  which  indicated  a  0.0048  atmospheres  (3.7  Torr)  of  water  vapor  and  an  ambient 
temperature  of  7°C.  As  evident,  the  DIAL  deduced  value  was  a  factor  of  2  less  than  the  actual 
experimental  conditions  indicated.  This  discrepancy  is  attributable  to  the  finite  linewidth  of  the 
Co:MgF2  laser  (0.15  cm"1)  in  comparison  to  the  water  vapor  linewidth,  which  serves  to  reduce 
the  peak  absorption  strength  and  to  broaden  the  measured  absorption  line,  as  seen  in  Figure  2. 
Cahen  and  Megie11  have  considered  this  problem  assuming  Gaussian  lineshapes.  They  found  that 
for  laser  linewidths  approximately  equal  to  the  absorption  linewidth,  deviations  in  the  observed 
absorption  strength  of  the  order  of  30%  may  be  expected,  and  that  the  deviation  increases  with 
increasing  laser  linewidth  and  integrated  optical  density.  It  is  clear  that  under  experimental  condi¬ 
tions  in  which  these  effects  are  significant,  accurate  molecular  concentration  determination  can 
only  occur  when  such  an  effect  has  been  calibrated  or  when  the  laser  linewidth  is  significantly 
narrower  than  the  molecular  absorption  linewidth. 

Frequency  scans  have  also  been  made  over  spectral  regions  containing  absorption  lines  of  the 
2.0  band  of  HC1  and  the  2v 3  band  of  CH4.  During  the  scan,  simultaneous  measurements  were 
made  of  the  return  signal  from  the  hillside  3  km  distant  and  the  transmittance  of  the  laser  beam 
passing  through  an  absorption  cell  1-m  in  length  containing  a  known  concentration  of  either  HC1 
or  CH4.  The  results  of  such  a  simultaneous  measurement  with  15  Torr  of  HC1  in  air  at  atmos¬ 
pheric  pressure  in  the  absorption  cell  are  shown  in  Figure  3.  The  frequency  scan  in  Figure  3  cov¬ 
ers  2.3  cm-1  and  encompasses  the  1.7525  R(l)  line  of  the  2-0  band  of  HC1.  The  upper  section 
of  Figure  3  corresponds  to  an  R(l)  absorption  coefficient  of  0.43  (atm-cm)'1.  The  remote  sensing 
data  obtained  simultaneously  and  shown  in  the  lower  portion  of  Figure  3  shows  the  presence  of 
nearby  water  vapor  lines,  but  within  experimental  uncertainty  there  is  little  indication  HC1 
(<0.2  ppm)  in  the  atmosphere  between  the  hillside  and  the  laboratory;  the  experimental  accuracy 
was  primarily  limited  by  uncertainty  in  the  background  signal  level  shown  in  Figure  3(b),  which 
led  to  an  error  of  approximately  5%. 

A  similar  measurement  was  carried  out  over  the  spectral  region  which  included  the  1.6713  jum 
P(2)  absorption  line  of  the  2v 3  band  of  CH4.  The  laboratory  measurement  indicated  an  absorp¬ 
tion  coefficient  of  0.036  (cm-atm)-1  for  this  line,  and  no  (<2  ppm)  CH4  was  detected  in  the 
atmosphere. 
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IV.  RANGE  RESOLVED  AEROSOL  DIAL  MEASUREMENTS 


The  Co:MgF2  DIAL  system  has  sufficient  transmitted  laser  energy  and  optical  detection  sen¬ 
sitivity  to  measure  the  lidar  backscatter  from  natural  oecurring  aerosols  in  the  atmosphere,  and 
thereby  provides  the  capability  of  detecting  the  range-resolved  concentration  of  various  species  in 
the  atmosphere.  As  an  example,  Figure  4  shows  the  lidar  signal,  averaged  over  50  pulses,  as  a 
function  of  range  for  atmospheric  backscatter.  As  seen,  the  Co:MgF2  DIAL  system  is  capable  of 
detecting  aerosol  backscatter  at  ranges  out  to  3.8  km,  with  significant  signal  levels  at  ranges  of 
0.2  to  3  km. 

Preliminary  experiments  using  the  range-resolved  Co:MgF2  DIAL  system  were  conducted  in 
order  to  assess  its  capabilities  for  atmospheric  studies.  Figure  5  shows  a  25  pulse  average  of  the 
range  resolved  lidar  returns  as  the  laser  was  tuned  on-  and  off-resonance  through  a  strong  water 
vapor  line  near  1.7515  ^im  with  an  expected  absorption  coefficient  of  8.3  X  10~3  (atm-cm)"1;  this 
strong  absorption  line  of  H20  was  chosen  to  enhance  the  differential  absorption  at  close  ranges. 
One  deduces  from  the  observed  differential  signal  near  300  m  in  Figure  5  that  approximately 
1.7  Torr  of  water  vapor  was  present.  I  his  value  differs  significantly  from  the  2.5  Torr  value  mea¬ 
sured  using  a  wet-bulb/ dry-bulb  thermometer.  The  difference  is  similar  to  that  observed  in  the 
path-averaged  measurements  (Figure  2)  and,  as  in  that  case,  the  deviation  is  presumably  due  to 
the  finite  linewidth  of  the  laser.  Calibration  errors  in  the  assumed  absorption  strengths  of  the 
water  vapor  lines  and  background  continuum  may  also  be  contributing  to  this  effect. 

The  Co:MgF2  DIAL  system  was  also  used  to  detect  HC1  in  the  atmosphere  under  range- 
resolved  conditions.  Since  the  background  concentration  of  HC1  in  the  atmosphere  is  too  small 
to  be  measurable,  artificial  means  of  injecting  a  plume  of  HC1  into  the  atmosphere  was  used.  In 
order  to  reduce  the  total  amount  of  gaseous  HC1  injected  into  the  atmosphere,  a  semi-enclosed 
cell  was  used.  The  1-m  long,  0.8-m  diameter,  window-less  cell  was  placed  at  a  distance  of  200  m 
and  the  0.3-m  diameter  lidar  beam  passed  unobstructed  through  the  cell.  Figure  6  shows  the 
normalized  lidar  signal  obtained  after  HC1  was  injected  into  the  cell.  The  difference  between  the 
on-line  and  off-line  signals  indicates  that  approximately  40  ppm  of  HC1  was  present  over  the 
200-m  range  resolution  of  the  system.  The  range  resolution  is  determined  by  the  laser  pulse 
length  and,  as  seen  in  the  figure,  is  responsible  for  the  considerable  degree  of  spatial  averaging 
seen  in  the  data.  A  shorter  pulse  length  of  the  Co:MgF2  laser  would  permit  a  shorter  range  reso¬ 
lution  to  be  obtained. 

V.  DIAL  MEASUREMENT  SENSITIVITY 

It  is  of  interest  to  calculate  the  theoretical  measurement  sensitivity  of  the  Co:MgF2  DIAL 
system.  The  sensitivity  of  a  differential  absorption  remote  sensing  system  is  determined  by  the 
minimum  change  in  the  signal  return  which  can  be  detected.  The  ultimate  sensitivity  at  far  ranges 
occurs  when  this  minimum  observable  change,  AP,  is  equal  to  the  noise  equivalent  power  (NEP) 
of  the  detection  system.  In  that  case,  for  measurements  of  signal  returns  from  a  target  located  at 
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a  range  R,  the  minimum  average  detectable  average  concentration  over  the  distance  to  the  target 
is  given,  for  small  absorption,  by12 

_  (NEP)  77 R 

Nmin_  2paaAKPt  exp  (~2aR)  ’ 

where  p  is  the  target  reflectivity,  aa  is  the  absorption  coefficient  of  the  molecule  being  investi¬ 
gated,  A  is  the  telescope  area,  K  is  the  optical  efficiency  of  the  system,  Pt  is  the  transmitted 
power,  and  a  is  the  atmospheric  absorption  exclusive  of  the  investigated  molecule.  For  the  case 
of  range-resolved  aerosol  backscatter,  aa  is  replaced  with  aaAR/R  and  the  target  reflectivity,  p,  is 
replaced  with  ft AR,  where  AR  is  the  range  resolution  of  the  lidar  and  ft  is  the  aerosol  backscat¬ 
ter  coefficient. 


At  close  in  ranges,  Nmjn  is  further  limited  by  fluctuations  in  the  return  signals  to  the  min¬ 
imum  fractional  change  of  those  signals,  AP/Pr,  which  can  be  detected;  Pr  is  the  returned  lidar 
power.  In  that  case,  for  hard  target  backscatter,12 


Nr 


(AP/Pr) 

2aaR 


(2) 


To  estimate  the  sensitivity  of  the  present  Co:MgF2  system  for  the  column  content  measure¬ 
ment  of  HC1  on  the  basis  of  Eqs.  (1)  and  (2),  the  parameters  used  were:  NEP  =  1.5  X  10~9  W, 
p  =  0. 1,  aa  =  0.43  (cm-atm)"1,  K  =  0.1,  A  =  0.125  m2,  Pt  =  3  X  104  W,  and  a  =  0.03  km'1.  The 
resulting  values  of  Nmin  as  a  function  of  range  are  shown  in  Figure  7,  where  it  is  assumed  that 
signal  averaging  will  reduce  the  minimum  measurable  AP/Pr  to  2%  and  will  produce  a  5-fold 
improvement  in  sensitivity  in  the  (NEP)-limited  region.  Improvement  of  this  magnitude  is  well 
within  the  capability  of  signal  averaging  for  both  path-averaged12  and  range-resoloved  measure¬ 
ments.14  The  number  of  pulses  required  to  achieve  a  five-fold  reduction  of  the  standard  deviation 
will  vary  with  atmospheric  conditions,  pulse-repetition  frequency  and  wavelength,  but  would 
normally  be  expected  to  occur  after  averaging  over  the  order  of  50  to  100  pulses.  The  figure 
indicates  that  the  system  should  be  capable  of  detecting  path  averaged  concentrations  of  HC1 
between  0.1  and  1  part  per  million  over  ranges  to  about  10  km  using  a  10  mJ  pulse  energy. 


Figure  8  shows  the  results  of  a  similar  analysis  carried  out  for  range-resolved  aerosol 
measurements;  in  this  case  /?  was  estimated  from  the  relationship  ji  -  (0.017/  V)  m"1  sr1,15  where 
V  is  the  visibility  in  meters,  and  assumed  to  be  10  km.  As  seen  in  Figure  8,  the  results  indicate 
that  a  concentration  sensitivity  of  a  few  parts  per  million  to  a  range  of  approximately  2  km  may 
reasonably  be  expected. 


A  calculation  of  the  sensitivity  of  the  Co:MgF2  DIAL  system  for  the  remote  sensing  of  CH4 
is  similar  to  that  shown  above,  except  that  the  concentration  detection  sensitivity  of  CH4  is 
approximately  an  order  of  magnitude  smaller  than  that  shown  for  HC1  in  Figures  7  and  8  due 
to  a  correspondingly  smaller  absorption  coefficient,  aa,  for  CH4.  In  this  case,  one  should  be  able 
to  measure  CH4  concentrations  of  10  to  100  parts  per  million,  which  could  be  useful  for 
detecting  and  measuring  gas  leaks. 
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VI.  CONCLUSIONS 


The  Co:MgF?  laser,  which  is  continuously  tunable  between  L5  and  2.3  /i m,  has  been  shown 
to  be  a  useful  coherent  radiation  source  for  the  remote  sensing  of  atmospheric  constituents  with 
absorption  lines  in  this  spectral  region.  The  Co:MgF2  remote  sensing  system  has  functioned  both 
in  a  mode  in  which  it  was  continuously  tuned  through  absorption  lines  of  H20,  HC1  and  CH4, 
and  in  a  differential  absorption  mode  operating  on  and  off  the  resonance  frequency  of  an  H20 
absorption  line.  While  the  remote  sensing  experiments  did  not  observe  any  atmospheric  HC1  or 
CH4,  calculations  indicate  that  the  present  system  should  be  capable  of  detecting  atmospheric 
HC1  and  CH4  at  concentrations  of  the  order  of  1  and  10  parts  per  million,  respectively. 
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FIGURE  CAPTIONS 


Fig.  1.  Schematic  of  Co:MgF2  DIAL  system. 

Fig.  2.  Lidar  returns  backscattered  from  a  hillside  as  a  function  of  the  Co:MgF2  laser 
wavelength;  The  absorption  features  are  due  to  water  vapor  in  the  atmosphere. 

Fig.  3.  Wavelength  scan  of  Co:MgF2  DIAL  system  showing  both  lidar  returns  and 
laboratory  absorption  cell  signals  as  the  laser  wavelength  was  scanned  over  an 
absorption  line  of  HC1.  No  HC1  was  measured  in  the  atmosphere. 

Fig.  4.  Range-resolved  Co:MgF2  lidar  returns  from  atmospheric  aerosols  as  a  function 
of  range. 

Fig.  5.  Atmospheric  range-resolved  lidar  returns  as  the  Co:MgF2  laser  was  tuned  on- 
and  off-resonance  to  the  absorption  line  of  water  vapor  near  1.7515  jum. 

Fig.  6.  Atmospheric  range-resolved  DIAL  measurements  of  HC1  injected  into  a  semi- 
enclosed  cell  located  at  a  range  of  200  m 

Fig.  7.  Calculated  path-averaged  sensitivity  of  the  Co:MgF2  DIAL  system  for  the 
remote  sensing  of  atmospheric  HC1. 

Fig.  8.  Calculated  range-resolved  sensitivity  of  the  Co:MgF2  DIAL  system  for  the 
remote  sensing  of  atmospheric  HCL 
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Figure  1. 
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LIDAR  SIGNAL  RETURNS  (Arbitrary  Units) 


(a)  HILLSIDE  AT  RANGE  =  3  km 


(b)  HILLSIDE  AT  RANGE  =  6.7  km 
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Figure  2. 


22 


LIDAR  SIGNAL  TRANSMITTANCE 

(Arbitrary  Units)  (Arbitrary  Units) 


(a)  SIGNAL  THROUGH  ABSORPTION  CELL 


(b)  RETURN  FROM  HILLSIDE:  R  =  3  Km 
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Figure  3. 
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APPENDIX  B 


The  following  is  a  reprint  of  a  journal  article  published  in  IEEE  Journal  of  Quantum  Electronics, 
March  1985,  entitled  “Frequency  Tuning  Characteristics  of  a  Q-Switched  Co:MgF2  Laser.” 
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Frequency  Tuning  Characteristics  of  a  Q-Switched 

Co:MgF2  Laser 

S.  L0VOLD,  P.  F.  MOULTON,  D.  K.  K1LL1NGER,  and  N.  MENYUK 


Abstract -A  tunable  ^-switched  Co;MgF2  laser  has  been  developed 
for  atmospheric  remote  sensing  applications.  Frequency  tuning  is  pro¬ 
vided  by  a  quartz  etalon  and  a  specially  designed  three-element  biie- 
fringent  filter  covering  the  whole  gain  bandwidth  of  the  Co:MgF2  laser. 
The  laser  has  good  temporal  and  speetral  characteristics,  with  an  emis¬ 
sion  linewidth  of  approximately  3  GHz  (0.1  cm"1). 

1.  Introduction 

HE  strong  coupling  between  electronic  energy  levels  and 
the  lattice  in  various  transition-metal  doped  crystals  has 
permitted  the  development  of  lasers  tunable  over  broad  fre¬ 
quency  regions  in  the  near  infrared.  The  first  such  laser  to  be 
demonstrated  was  the  ^fi-doped  MgF3  laser  by  Johnson  etal. 
in  1963  [1].  Later,  laser  pumping  of  Ni-,  Co-  and  V-doped 
MgF2  crystals  greatly  improved  the  performance  of  these 
so-called  “vibronic”  lasers  [2] -[6] .  The  Co:MgF2  laser  is  one 
of  the  most  promising  among  these  because  it  has  negligible 
excited  state  absorption  and  can  be  tuned  from  1 .6  to  2.3  iw\ . 
CW,  ^-switched,  and  mode-locked  operation  has  been  demon¬ 
strated,  with  up  to  80  percent  quantum  efficiency  [3] -[6] . 

The  gain  curve  of  the  Co.MgF2  laser  has  several  sharp  peaks 
at  pure  electronic  transitions,  as  seen  in  Fig.  1.  The  peak  gain 
cross  section  is  relatively  low,  1.5  X  10_2lcm2  [7].  Thus  the 
precise  tuning  of  the  laser  across  the  gain  bandwidth  requires  a 
low  loss  filter  with  a  narrow  bandwidth  to  avoid  frequency 
pulling  of  the  laser  frequency  towards  the  peaks  in  the  gain 
curve. 

A  birefringent  filter  has  the  possibility  of  meeting  the  above 
requirements.  The  theory  and  applications  of  birefringent  fil¬ 
ters  have  been  extensively  covered  in  the  literature  [8] -[14] . 
In  order  to  provide  unambiguous  continuous  frequency  tuning 
across  the  whole  gain  bandwidth  of  the  Co:MgF2  laser, 
without  possible  frequency  jumps  between  various  filter  res¬ 
onances,  one  should  choose  the  optical  axis  of  the  bire¬ 
fringent  plates  to  make  an  angle  to  the  surface  of  the  plates. 
This  case  has  not  been  covered  in  detail  in  the  literature. 

In  this  paper  an  analysis  of  a  multielement  birefringent  filter 

Manuscript  received  July  9, .1984,  revised  November  12,  1984  This 
work  was  supported  in  part  by  the  National  Aeronautical  and  Space 
Administration,  by  the  Air  Force  Engineering  and  Service  Center,  and 
by  the  Air  Force  Aerospace  Medical  Division.  S.  L^vold  was  also  sup¬ 
ported  in  part  by  the  Norwegian  Defence  Research  Establishment  and 
by  the  Royal  Norwegian  Council  for  Industrial  and  Scientific  Research. 

S.  L/vold  was  with  the  Lincoln  Laboratory,  Massachusetts  Institute 
of  Technology,  Lexington,  MA  02173.  He  is  now  with  the  Norwegian 
Defence  Research  Establishment,  N-2007,  Kjeller,  Norway. 

P.  F.  Moulton,  D.  Killinger,  and  N.  Menyuk  are  with  the  Lincoln 
Laboratory,  Massachusetts  Institute  of  Technology,  Lexington,  MA 
02173. 


WAVELENGTH  (pm) 

Fig.  1.  Measured  absorption  and  fluorescence  spectrum  for  Co:MgF2 
between  1.0  and  2.3  um. 

is  presented  which  includes  the  case  where  the  optical  axis 
makes  an  angle  to  the  plate  surface.  This  analysis  constitutes 
a  major  part  of  the  paper.  A  three-element  filter  has  been  de¬ 
signed  [15]  which  satisfies  the  requirements  listed  above  for 
the  Co  MgF2  laser.  The  first  high-frequency-resolution  tuning 
characteristics  of  a  ^-switched  Co:MgF2  laser  are  presented. 
This  laser  has  been  developed  for  atmospheric  spectroscopy 
and  remote  sensing  applications,  which  require  that  the  laser 
linewidth  be  less  than  or  equal  to  3  GHz  (0.1  cm"1),  corre¬ 
sponding  to  the  pressure  broadened  linewidth  at  1  atm  of  most 
target  species  of  interest.  This  goal  is  achieved  by  insertion  of 
an  etalon  into  the  laser  cavity  in  addition  to  the  birefringent 
filter. 

II.  Theory 

A  birefringent  filter  for  laser  frequency  tuning  consists  of 
one  or  more  birefringent  plates  oriented  at  Brewster’s  angle  in¬ 
side  the  laser  resonator.  The  plates  act  both  as  retarding  plates 
and  polarizing  elements,  p-polarized  light  incident  on  the 
plates  is  transformed  into  ehptically  polarized  light  inside  the 
plates  with  continuously  varying  eccentricity  of  the  ellipse  as 
the  beam  propagates  through  the  crystal.  Certain  wavelengths 
will  be  linearly  p-polarized  at  the  exit  surface  of  the  plates  and 
will  be  transmitted  with  no  Fresnel  reflection  losses.  The  filter 
is  tuned  by  rotating  the  plates  about  the  axis  normal  to  the 
surface,  thereby  varying  the  effective  principal  refractive  in¬ 
dexes  of  the  plates. 

The  main  reason  for  choosing  the  optical  axis  to  make  an 
angle  to  the  plate  surface  is  to  increase  the  tuning  rate  versus 
plate  rotation.  This  is  important  because  of  the  narrowness  of 
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Fig.  2.  A  one-clement  birefringent  filter.  The  plate  is  oriented  at 
Brewster’s  angle  to  the  incoming  beam.  0  Is  the  internal  Brewster’s 
angle  and  s  is  the  direction  of  beam  propagation  in  the  plate.  The  c- 
axis  is  the  optical  axis  of  the  plate.  The  unit  vectors  ep  and  es  are 
parallel  to  the  p-  and  j-polarization  of  light  incident  on  the  plate.  The 
plate  is  rotated  an  angle  p  defined  in  the  figure,  y  is  the  angle  be¬ 
tween  the  c-axis  and  the  direction  of  beam  propagation,  and  r  is  the 
thickness  of  the  plate. 

the  plate  rotation  domains  where  the  modulation  depth  of  the 
Filter  is  good. 

In  addition,  one  of  these  domains  is  broader  when  the  opti¬ 
cal  axis  points  out  of  the  surface,  compared  to  when  it  lies  in 
the  plane  of  the  plate.  These  effects  will  be  shown  in  the  fol¬ 
lowing  discussion. 

A  one-element  birefringent  Filter  is  shown  schematically  in 
Fig.  2.  The  light  is  incident  on  the  plate  at  Brewster’s  angle. 
(0  is  the  internal  Brewster’s  angle.)  Using  the  Jones  vector  for¬ 
malism  (16),  the  electric  field  vector  £  is  represented  as 

E  -  Epep  +  Ej,  (1) 

where  Ep  and  Ea  are  the  p-  and  s-polarized  electric  field  vector 
components  along  the  unit  vectors  ep  and  e3  respectively  The 
c-axis  of  the  plate  (the  optical  axis)  makes  an  angle  a  to  the 
surface  normal,  and  the  plate  is  rotated  an  angle  p  as  defined 
in  the  Figure. 

After  transmission  through  the  plate,  the  electric  field  com¬ 
ponents  are  given  by 


where  the  matrix  I  describes  the  transmission  through  a  plate 
interface  and  the  matrix  M  describes  the  transmission  through 
the  birefringent  material.  (•  means  standard  matrix  multiplica¬ 
tion.)  1  is  given  by  (17) 


where  q  -  7nl(n2  +  1)  is  the  transmittance  of  the  s-polarized 
electric  field  at  Brewster’s  angle  into  a  medium  with  index  of 
refractions.  The  matrix  M  is  given  by  [16],  [11] 

/cos2  0  +  sin2  0e/A<?>  cos  0  sin  0(e;A<?>  -  1)\ 

M  =  l  (4) 

\cos  0  sin  d(e,A<t>  -  1)  sin2  0  +  cos2  0e/  A<?  / 


Fig.  3.  The  index  ellipsoid  of  a  positive  uniaxial  birefringent  crystal 
[171.  The  vector  s  is  parallel  to  the  direction  of  the  beam  propa¬ 
gation  in  the  crystal,  which  makes  an^angle  7  to  the  optical  axis  (the 
c-axis).  The  orthogonal  unit  vectors  ep  and  e3  are  parallel  to  the  p- 
and  j-polanzation  of  the  electric  field,  inside  the  plate.  The  ortho¬ 
gonal  unit  vectors  eQ  and  ee  are  parallel  to  the  semiminor  and  semi¬ 
major  axis,  respectively,  of  the  ellipse  defined  by  the  intersection  of 
the  j -plane  with  the  ellipsoid,  d  is  the  angle  between  ep  and  eQ. 

Here  0  is  the  angle  between  the  unit  vectors  e0  and  Vp.  The 
vector  eQ  is  parallel  to  the  direction  in  the  plate  where  the 
index  of  refraction  is  n0,  and  ep  and  es  are  the  unit  vectors 
parallel  to  the  p-  and  s -polarization  for  the  electric  field  inside 
the  plate,  as  shown  in  Figs.  2  and  3.  (The  matrix  M  is  ob¬ 
tained  by  simply  transforming  the  electric  field  from  the 
coordinate  system  ep,  e'3  to  e0l  ee  (see  Fig.  3),  and  back 
again  to  ep,  P3  after  propagation  through  the  crystal.)  A0  is 
the  phase  retardation  of  the  plate,  and  is  given  to  a  good  ap¬ 
proximation  when  ne  -  nQ  «  nQ  by  [17] 

M  =  («»  -  n0)  sin2  r  (5) 

A  COS  p 

Here  t  is  the  thickness  of  the  plate,  nQ  and  n€  are  the  ordinary 
and  extraordinary  indexes  of  refraction  respectively,  and  y  is 
the  angle  between  the  c-axis  and  the  direction  of  the  beam 
propagation  5  in  the  crystal.  The  wavelengths  for  which  A <p  = 
N2n ,  where  N  is  an  integer,  will  have  the  same  polarization  at 
the  exit  interface  as  at  the  entrance  interface  of  the  plate. 

The  angles  7  and  0  entering  (4)  and  (5)  can  be  expressed  by 
the  angles  p,  and  o  (Fig.  2)  through  the  equations 

cos  7  =  cos  (3  cos  o  +  sin  (3  sin  o  cos  p  (6) 

cos  0  =  sin  p  sin  0/sin  7.  (7) 

(Eqs.  ^6)  and  (7)  can  be  obtained  from  the  vector  multiplica¬ 
tions  c  •  5  =  ici  •  IjI  cos  7  and  ep  •  s  X  c  =  Isl  •  lei  sin  7  cos  0, 
as  can  be  seen  from  Figs.  2  and  3.) 

With  the  possible  addition  of  an  equation  describing  the  dis¬ 
persion  of  the  birefringent  material,  (l)-(7)  form  the  basis  for 
designing  a  tuning  filter  for  a  given  laser.  For  crystal  quartz, 
the  dispersion  equation 

ne-  n0=  [8.5  -  0.65(X  -  1 .5)]  X  1CT3  (8) 

is  a  good  approximation  to  experimental  data  for  wavelengths 
X  (in  micrometers)  between  0  6  and  2  .4  pm  [18]. 

In  designing  a  Filter  for  a  given  wavelength  region,  care 
should  be  taken  that  the  resonance  condition,  A <p  =  N2n,  does 
not  occur  for  plate  rotation  angles  where  the  electric  field 
polarization  is  parallel  to  the  ordinary  or  extraordinary  axis  of 
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Fig.  4.  The  modulation  depth,  defined  by  (9),  versus  plate  rotation 
angle  p  for  various  c-axis  orientations,  a  is  the  angle  between  the 
c-axis  and  the  plate  normal  (see  Fig.  2). 


Fig.  5.  Resonance  wavelengths,  lor  which  the  transmission  for  incident 
p-polarized  light  is  100  percent,  versus  plate  rotation  angle,  for  two 
different  single-element  filters.  One  filter  is  a  0.65  mm  thick  plate 
with  angle  a  *  35°  between  the  optical  axis  and  the  plate  normal. 
The  resonances  are  indicated  by  solid  lines.  The  other  filter  is  a  0.68 
mm  thick  plate  with  o  =  90°,  i.e.,  the  optical  axis  is  in  the  plane  of 
the  plate.  The  resonances  are  indicated  by  the  dashed  curves.  Both 
plates  arc  of  single  crystal  quartz. 

the  crystal.  In  these  cases  all  wavelengths  will  be  transmitted 
with  no  Fresnel  reflection  losses.  This  corresponds  to  the 
situation  where  the  angle  d  (Fig.  3)  is  0  or  90°.  It  is  useful 
to  define  the  “modulation  depth1’  MD  by 

MD  =  \m2l  (max)l2  =  4  cos2  d  sin2  d  (9) 

where  m2 1  is  the  off-diagonal  element  of  the  matrix  M  {m2x  - 
m12).  Fig.  4  shows  the  modulation  depth  versus  plate  rotation 
angle  for  various  c-axis  orientations.  The  domains  where  MD 
is  large  defines  the  useful  plate  rotation  angle  domains.  A 
design  criterion  for  the  filter  is  thus  that  only  one  filter  reso¬ 
nance  falls  within  such  a  domain  and  that  this  resonance 
covers  the  whole  gain  bandwidth  of  the  laser  in  that  domain. 

Fig.  5  shows  the  resonance  wavelengths  for  two  different 
plates.  One  is  a  0.65  mm  thick  crystal  quartz  plate  with  an 
angle  a  =  35°  between  the  c-axis  and  the  surface  normal.  The 
other  is  0.68  mm  thick  of  the  same  material  with  a  =  90°,  i.e., 
the  c-axis  lies  in  the  plane  of  the  plate.  It  is  clearly  seen  that 
the  plate  with  a  =  35°  can  be  useful  for  tuning  a  Co:MgF2 
laser  from  1 .5  to  2.3  pm  on  the  N  s  2  resonance.  This  reso¬ 
nance  falls  within  a  domain  where  the  modulation  depth  is 
large,  as  seen  in  Fig.  4.  At  the  longer  wavelengths  the  inter¬ 
ference  from  the  resonance  N  -  3  can  cause  a  problem.  How¬ 
ever,  the  laser  mirrors  cannot  easily  be  made  to  cover  the 
whole  spectral  range  from  1.5  to  2.3  pm  with  uniform  reflec¬ 
tivity.  Thus  this  possible  interference  should  not  cause  any 
problems.  If  there  should  be  an  interference  problem  between 
the  A^  =  2  and  N  =  3  resonances,  one  could  use  an  additional 
plate  with  half  the  thickness,  i.e.,  (  =  0,325  mm.  This  plate 
would  have  its  N  ~  1  resonance  overlapping  the  N  =  2  reso¬ 


nance  of  the  0.65  mm  plate,  and  would  have  large  losses  at  the 
N  -  1  and  N  88  3  resonances  of  the  0.65  mm  plate. 

It  can  be  seen  in  Fig.  4  that  when  the  optical  axis  lies  in  the 
plane  of  the  plate,  i.e.,  o  =  90°,  a  single  resonance  order  is 
limited  to  the  spectral  region  from  1.6  to  2.15  pm  when  the 
plate  is  rotated  from  0  to  90°.  In  addition,  it  was  shown  in 
Fig.  4  that  the  modulation  depth  is  poor  close  to  0  and  90°. 
This  is  therefore  not  a  good  solution.  It  should  be  noted  that 
one  could  choose  a  thicker  plate  and  use  two  or  more  reso¬ 
nance  orders  to  cover  the  gain  bandwidth  of  the  Co:MgF2 
laser.  This  would  require  a  set  of  mirrors  with  center  fre¬ 
quencies  and  bandwidths  that  were  matched  to  such  a  filter. 
In  this  work,  however,  the  filter  with  o  =  35°  has  been  found 
to  be  a  better  solution. 

The  transmission  bandwidth  of  the  filter  (for  a  fixed  rota¬ 
tion  angle)  determines  the  degree  of  laser  frequency  pulling 
towards  peaks  in  the  gain  curve,  and  also  determines  the  laser 
emission  linewidth. 

The  transmission  T  of  a  birefringent  filter  can  be  defined  as 
_  Ep,  out 

T=  - 

E  • 

up ,  in 

that  is,  the  fraction  of  the  incident  p-polarization  intensity 
that  is  left  in  the  ^-polarization  after  propagation  through 
the  filter.  If  the  incident  light  is  fully  p-polarized  and  the  bi¬ 
refringent  plates  were  followed  by  a  100  percent  p-polarizer,  T 
would  be  the  actual  transmission  of  the  filter.  We  assume  that 
the  transmission  T  as  defined  by  (10)  is  useful  for  characteriz¬ 
ing  the  transmission  of  a  birefringent  filter  for  a  Co:MgF2 
laser  because  the  modes  with  the  highest  net  gam  in  such  a 
laser  will  be  linearly  polarized.  This  is  due  to  the  birefrin¬ 
gence  of  the  MgF2  crystal  and  usually  the  presence  of  several 
polarizing  Brewster’s  angle  interfaces  in  the  resonator  (i.e.,  the 
MgF2  crystal  interfaces,  the  ^-switch,  and  the  vacuum-Dewar 
windows-  see  Section  111). 

The  transmission  bandwidth  of  a  single  plate  filter  operating 
on  a  low-order  resonance  will  be  broad,  typically  several  hun¬ 
dred  nanometers.  Adding  more  plates  with  c-axes  parallel  to 
each  other  and  thickness  equal  to  a  multiple  of  the  thickness 
of  the  first  plate  will  make  the  transmission  bandwidth  nar¬ 
rower,  but  will  also  introduce  secondary  transmission  peaks. 
Fig.  6  shows  the  transmission  T  as  defined  by  (9)  fora  three- 
element  birefringent  filter  with  c-axis  angle  o  =  35°  and  thick¬ 
ness  /,  2 1,  and  1  Or ,  where  t  =  0.65  mm.  The  plate  rotation 
angle  p  =  90°.  The  peak  transmission  wavelength  is  at  the 
second  order  resonance  in  the  thinnest  plate  and  at  the  20th 
order  resonance  in  the  thickest  plate. 

In  order  to  evaluate  whether  the  laser  wavelength  will  follow 
the  main  transmission  peak  of  the  filter,  the  transmission  loss 
on  the  secondary  peaks  should  be  compared  with  the  maxi¬ 
mum  difference  in  single-pass  gain  in  the  gain  medium  over 
corresponding  wavelength  intervals.  For  the  filter  shown  in 
Fig.  6  the  transmission  losses  are  15  percent  at  the  secondary 
transmission  peaks,  at  the  19th  and  21st  order  resonance  of 
the  thickest  plate  For  the  Co:MgF2  laser  the  peak-to-valley 
gain  ratio  over  a  corresponding  wavelength  interval  is  7/5  at 
around  1 .7  pm,  according  to  the  fluorescence  curve  shown  in 
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Fig.  6.  The  transmission  T,  as  defined  by  (10)  for  a  three-element 

crystal  quartz  birefringent  filter.  The  c-axis  angle  o  -  35°.  The  c- 

axes  are  parallel  to  each  other  and  the  thicknesses  of  the  plates  are  t. 

2 1,  and  10?,  where  t  a  0.65  mm.  The  plate  rotation  angle  p  =  90°. 

Fig.  1 .  Assuming  a  laser  resonator  with  an  average  4  5  percent 
single-pass  loss  (for  example  7  percent  output  coupling,  1  per¬ 
cent  single-pass  internal  losses),  one  can  calculate  that  the  laser 
can  be  pumped  3  times  above  threshold  for  operation  at  a 
wavelength  in  the  gain  valley  domain  before  the  net  laser  gain 
can  be  positive  at  the  peak  gain  wavelength.  The  laser  can  be 
pumped  8  times  above  threshold  for  operation  in  the  gain 
valley  domain  before  the  net  gain  at  the  peak  can  exceed  the 
net  gain  in  the  valley  .  The  filter  shown  in  Fig.  6  has  more  than 
13  percent  loss  at  the  secondary  transmission  peaks  across  the 
whole  gain  bandwidth  of  the  Co:MgF2  laser.  Thus  since  the 
Co:MgF2  is  such  a  low  gain  material  and  it  is  predominantly 
homogeneously  broadened  [4] ,  this  filter  should  have  an  ade¬ 
quate  extinction  ratio. 

The  laser  emission  linewidth  is  determined  by  the  transmis¬ 
sion  bandwidth  of  the  filter,  the  mode  of  operation  of  the 
laser,  and  whether  the  laser  transition  is  homogeneously  or 
inhomogeneously  broadened.  In  ^-switched  operation  the 
linewidth  will  be  determined  by  the  filter  transmittance  raised 
to  the  power  of  2A'rt,  where  Nn  is  the  number  of  round  trips 
a  photon  makes  in  the  resonator  during  the  build-up  time  of 
the  laser  pulse  In  the  Co:MgF2  laser  used  in  our  experiments, 
iVrt  *  400  when  the  laser  is  pumped  1.5  times  above  thresh¬ 
old,  according  to  0-switched  laser  theory  {19],  [20] .  The 
transmittance  of  the  filter  shown  in  Fig.  6,  raised  to  the  power 
of  2Nn,  will  have  a  central  peak  width  (FWHM)  of  115  GHz 
(3.8  cm”1)  for  Nrt  =  400.  Since  the  switching  time  for  the 
acoustooptic  0-switch  used  in  these  experiments  is  relatively 
long  (200-500  ns,  corresponding  to  40-100  round-trips)  and 
since  the  laser  transition  is  predominantly  homogeneously 
broadened,  the  laser  emission  linewidth  is  expected  to  be  nar¬ 
rower,  but  of  the  same  order  of  magnitude,  as  the  calculated 
1 15  GHz,  for  pumping  1.5  times  above  threshold.  This  is  in 
good  agreement  with  the  experimental  results  (see  Section  IV). 

The  number  of  round  trips  during  the  build-up  time  of  the 
0-switched  laser  pulse  is  roughly  proportional  to  {Tq  [(P/Pt)  " 

1  ]  }“ 1  where  T0  is  the  output  coupling  and  P/Pt  is  the  pumping 
relative  to  the  threshold  pump  level  [19]  [20] .  Since  the  “ef¬ 
fective”  transmission  bandwidth  of  the  filter  is  inversely  pro- 


Fig.  7.  The  Co  MgF2  laser  system.  The  Co:MgF2  crystal  is  mounted 
inside  a  liquid  nitrogen  Dewar.  A  three-clement  birefringent  filter 
(BRF)  and  a  quartz  etalon  (ET)  are  used  for  tuning  of  the  Co  :MgF2 
laser,  and  the  laser  has  an  acoustooptic  0-switch  (AOQS).  The 
Co:MgF2  laser  is  pumped  by  a  Nd:YAG  laser  with  an  internal 
Brewster  window  (BW),  to  force  it  to  emit  linearly  polarized  light. 

portional  to  the  square  root  of  the  number  of  transmissions 
through  the  filter,  it  is  expected  that  the  laser  linewidth  will 
vary  accordingly  when  the  pump  level  or  the  output  coupling 
are  changed. 

111.  The  Co:MgF2  Laser 

The  Co:MgF2  laser  system  is  shown  schematically  in  Fig.  7. 
The  Brewster-angle-cut  MgF2  crystal  is  27  X  6  X  6  mm  in  size 
and  doped  with  0.75  weight  percent  CoF2.  It  is  mounted 
inside  a  liquid  nitrogen  Dewar  with  Brewster  windows  and 
cooled  to  77  K.  The  resonator  mirror  spacing  is  70  cm  and  the 
mirror  radii  arc  100  cm.  The  output  mirror  has  5  percent 
transmission,  centered  at  1.75  /im.  The  Co:MgF2  laser  is  lon¬ 
gitudinally  pumped  by  a  Nd:YAG  laser,  which  emits  400  mJ, 
800  {is  duration  linearly  polarized  multimode  pulses  at  1.32 
Mm.  (Without  the  Brewster  window  inside  the  Nd:  YAG  laser 
resonator  ,  the  output  energy  is  500  mJ  randomly  polarized 
for  the  same  input  energy.)  The  pump  pulses  are  focused 
using  a  25  cm  focal  length  lens  to  a  1  mm  diameter  spot  in  the 
Co  MgF2  laser  crystal,  exciting  only  the  TEM0o  mode  of  the 
Co :  MgF2  laser  resonator. 

The  laser  frequency  can  be  tuned  by  a  three-element  bire¬ 
fringent  filter  (BIFF  or  BRF)  and  a  1  mm  thick  uncoated  quartz 
etalon.  The  BRF  has  been  chosen  with  c-axis  angle  a  =  35° 
and  thickness  t,  2 r,  and  1  Or  where  t  -  0.65  mm,  as  described 
in  the  preceeding  section.  Alignment  of  the  individual  plates 
with  their  c-axis  parallel  to  each  other  was  first  done  external 
to  the  laser  using  a  Hc-Ne  laser  and  calculated  rotation  angles 
for  resonances  at  0.6328  Mm.  The  filter  was  then  mounted 
inside  the  Co  MgF2  laser  and  the  individual  plates  were  fine- 
adjusted  for  maximum  laser  output.  From  measurements  of 
the  threshold  for  laser  oscillation  before  and  after  insertion  of 
the  filter,  the  single-pass  insertion  loss  was  calculated  to  be  less 
than  1  percent 

The  laser  was  0-switched  by  a  Brewster-angle-cut  fused 
silica  acoustooptic  0-switch.  The  single-pass  insertion  loss  of 
the  0-switch  was  measured  to  be  less  than  0.5  percent.- 

IV.  Results 

The  temporal  and  spectral  characteristics  of  the  laser  output 
have  been  measured.  The  spectral  characteristics  were  mea¬ 
sured  with  a  1  m  grating  spectrometer  with  3  GHz  (0.1  cm-1 ) 
resolution,  and  with  a  plane-parallel  scanning  Fabry-Perot 
interferometer.  The  system  finesse  of  the  Fabry-Perot  inter¬ 
ferometer  has  been  measured  to  be  better  than  20.  The  laser 
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Fig.  8.  The  measured  laser  output  wavelength  versus  biicfringent 
filter  rotation,  and  the  measured  output  pulse  energy  versus  wave¬ 
length.  The  laser  was  operated  in  a  free-running  (gain-switched)  mode 
without  the  0-switch  and  the  etalon.  The  transmission  losses  of  the 
laser  mirrors  are  also  shown. 

frequency  has  also  been  scanned  across  an  absorption  line  of 
HC1. 

Fig.  8  shows  the  measured  laser  output  wavelength  versus 
BRF  rotation,  and  the  measured  output  pulse  energy  versus 
wavelength.  In  this  case  the  0-switch  and  the  etalon  were 
removed.  The  laser  was  pumped  with  approximately  400  mJ 
at  a  3  Hz  pulse  repetition  rate.  It  can  be  seen  that  the  laser 
wavelength  tunes  smoothly  with  the  BRF  rotation,  and  there 
is  an  excellent  agreement  between  the  measured  tuning  charac¬ 
teristics  and  those  calculated  for  N  =  2,  as  shown  in  Fig.  5. 
No  frequency  pulling  towards  the  peaks  in  the  gain  curve 
could  be  observed  to  within  the  measurement  accuracy  of  the 
grating  spectrometer.  The  same  smooth  tuning  characteristics 
were  obtained  during  0-switched  operation. 

With  the  present  set  of  resonator  mirrors,  the  laser  could  be 
tuned  from  1 .6  to  1 .9  iim.  At  wavelengths  longer  than  1 .8  jrm 
the  combined  effect  of  mirror  loss  and  strong  absorption  from 
water-vapor  in  the  air  limits  the  output  energy  and  tuning 
range.  The  total  transmission  of  the  two  resonator  mirrors  is 
shown  in  Fig.  8.  (The  water-vapor  absorption  can  of  course  be 
avoided  by  placing  the  whole  resonator  in  vacuum  or  in  a  dry 
nitrogen  atmosphere.  This  has  not  been  done  in  these  experi¬ 
ments.)  The  two  strong  peaks  in  the  output  energy  at  1 .67 
and  1.75  nm  correspond  to  electronic  transitions  in  the  gain 
spectrum,  as  seen  in  Fig.  1. 

The  temporal  and  spectral  characteristics  of  the  laser  output 
were  measured  during  0-switched  operation.  The  laser  was 
pumped  approximately  1.5  times  above  threshold  at  various 
wavelengths  between  1.67  and  1.8  jim.  The  measured  output 
pulse  energy  and  pulse  duration  was  ~15  mi  and  600  ns 
respectively,  in  good  agreement  with  0-switched  laser  theory 
[19], [20]. 

In  0-switched  operation  and  without  an  etalon  in  the  cavity 
the  laser  emission  linewidth  was  measured  to  be  ^100  GHz 
(~3  cm  “ 1 ).  This  is  in  good  agreement  with  what  could  be  ex¬ 
pected  from  the  transmittance  characteristics  of  the  BRF,  as 
discussed  in  Section  11. 

When  the  1  mm  thick  uncoated  quartz  etalon  was  inserted 


Fig.  9.  Fabry-Perot  interferometer  scan  of  the  laser  beam,  for  a  fixed 
laser  frequency.  The  free  spectral  range  of  the  interferometer  is  15 
GHz  (1  cm  mirror  separation).  The  laser  PRF  is  3  Hz  and  the  laser 
pulses  are  sampled  by  a  boxcar  signal  averager.  The  interferometer 
scan  rate  is  12  GHz/ min. 

into  the  cavity,  the  laser  emission  linewidth  as  measured  with 
the  scanning  Fabry-Perot  interferometer  was  found  to  be  3 
GHz  (0  1  cin“‘).  Fig  9  shows  an  interferometer  scan  with  a 
1 5  GHz  free  spectral  range  (1  cm  mirror  separation)  for  a  fixed 
laser  frequency.  The  frequency  scan  is  very  slow,  so  the 
spectrum  is  an  average  taken  over  many  laser  pulses.  Appar¬ 
ently  the  laser  oscillates  on  several  longitudinal  modes  simul¬ 
taneously.  The  pulse-to-pulse  frequency  variation  and  the  fre¬ 
quency  spectrum  of  a  single  pulse  have  not  been  measured. 

The  measured  3  GHz  laser  emission  linewidth  is  close  to  the 
value  expected  on  the  basis  of  the  transmission  characteristics 
of  the  etalon.  Assuming  400  round-trips  in  the  resonator 
during  the  build-up  time  for  the  laser  pulse  for  pumping  1 .5 
times  above  threshold,  the  bandwidth  of  the  etalon  transmis¬ 
sion  lobes  will  be  5  GHz  (FWHM)  (see  Section  11). 

A  laser  frequency  scan  requires  a  simultaneous  tuning  of  the 
birefringent  filter  and  the  etalon.  The  laser  may  be  tuned 
across  a  10  GHz  frequency  range  (the  free  spectral  range  of  the 
etalon)  by  tilting  the  etalon  and  keeping  the  BRF  fixed.  Fig 
10  shows  such  a  scan  across  the  first  overtone  of  the  R3  line  of 
HC137,  at  X  =  1 .7436  /urn  [21  ] .  The  absorption  cell  contained 
50  torr  of  HCl  and  700  torr  of  air.  The  measured  linewidth  is 
approximately  5  GHz  which  agrees  well  with  the  expected 
pressure  broadened  linewidth  at  1  atm  pressure  During  the 
scan  shown  in  Fig.  10  the  laser  was  operated  in  a  free-running 
(gain  switched)  mode,  in  which  case  the  laser  linewidth  has 
been  measured  to  be  less  than  1 .2  GHz  (0.04  cm  " 1 ). 

V.  Discussion  and  Conclusion 
In  order  to  keep  the  output  energy  level  well  below  the 
threshold  for  laser  mirror  damage  [6] ,  the  laser  was  pumped 
only  1.5  times  above  threshold  in  these  experiments.  A  lidar 
system  has  been  built,  using  this  laser  as  a  transmitter,  which 
has  demonstrated  path-averaged  absorption  spectroscopy  of 
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Fig.  10.  A  laser-frequency  scan  across  the  first  overtone  of  the  R3  line 
of  HC137,  at  171.9  TH z  (X  *  1.7436  nm).  The  laser  was  operated  in 
a  free-running  (gain-switched)  mode.  The  absorption  cell  contained 
50  torr  of  HC1  and  700  torr  of  air. 

atmospheric  species  with  3  GHz  spectral  resolution  out  to  7.35 
km  [22] .  Separate  experiments  have  demonstrated  that  by 
using  a  higher  output  coupling,  at  least  60  mJ  output  pulses 
with  150  ns  duration  can  be  extracted  from  the  Co:MgF2  laser 
at  1.96  pm  without  damage  to  the  laser  mirrors  or  the  crystal. 
This  should  allow  range-resolved  atmospheric  remote  sensing 
measurements  to  be  made  with  25  m  spatial  resolution  out  to 
several  kilometers  [23) ,  [24 ] . 

It  should  be  noted  that  when  the  laser  is  pumped  further 
above  threshold  or  when  the  output  coupling  is  increased,  the 
build-up  time  for  the  ^-switched  laser  pulse  is  decreased.  This 
is  expected  to  give  a  broader  laser  emission  linewidth,  as  dis¬ 
cussed  in  Section  11  It  may  thus  then  be  necessary  to  use  a 
sharper  filter  than  the  1  mm  thick  uncoated  quartz  etalon  in 
order  to  obtain  a  3  GHz  or  smaller  laser  emission  linewidth. 
As  an  alternative,  the  loss  of  the  (9-switch  can  be  adjusted  to 
allow  the  laser  to  operate  slightly  above  threshold  during  the 
initial  stage  of  the  pump  pulse,  thus  allowing  formation  of  a 
narrow-line  width  signal  in  the  cavity.  The  ^-switch  loss  can 
then  be  reduced  to  zero  to  generate  a  short,  energetic  pulse. 

In  conclusion,  a  tunable  ^-switched  Co:MgF2  laser  system 
has  been  developed  which  has  demonstrated  that  the 
Co:MgF2  laser  has  good  temporal  and  spectral  characteristics 
which  are  suitable  for  atmospheric  spectroscopy  and  remote 
sensing  applications. 
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Laser  Remote  Sensing  of  the  Atmosphere 

Dennis  K.  Killinger  and  Norman  Menyuk 


Laser  beams  can  be  used  as  long-range  spectroscopic 
probes  of  the  chemical  composition  and  physical  state  of 
the  atmosphere.  The  spectroscopic,  optical,  and  laser 
requirements  for  atmospheric  laser  remote  sensing  are 
reviewed,  and  the  sensitivity  and  limitations  of  the  tech¬ 
nique  are  described.  A  sampling  of  recent  measurements 
includes  the  detection  of  urban  air  pollution  and  toxic 
chemicals  in  the  atmosphere,  the  measurement  of  global 
circulation  of  volcanic  ash  in  the  upper  atmosphere,  and 
the  observation  of  wind  shear  near  airports. 


Laser  remote  sensing  of  atmospherjc  properties 
from  a  single  location  is  referred  to  as  lidar,  an  acronym  for 
light  detection  and  ranging,  and  is  analogous  to  radar  (1-6). 
In  lidar,  the  projection  of  a  short  laser  pulse  is  followed  by  reception 
of  a  portion  of  the  radiation  reflected  from  a  distant  target  or  from 
atmospheric  constituents  such  as  molecules,  aerosols,  clouds,  or  dust 
(Fig-  1)-  The  incident  laser  radiation  interacts  with  these  constitu¬ 
ents,  causing  alterations  in  the  intensity  and  wavelength  according 
to  the  strength  of  this  optical  interaction  and  the  concentration  of 
the  interacting  species  in  the  atmosphere.  Conscquendv,  informa¬ 
tion  on  the  composition  and  physical  state  of  the  atmosphere  can  be 
deduced  from  the  lidar  data.  In  addition,  the  range  to  the  interacting 
species  can  be  determined  from  the  temporal  delay  of  the  backseat* 
tered  radiation.  Lidar  has  been  used  to  measure  the  movement  and 
concentration  of  air  pollution  near  urban  centers,  the  chemical 
emission  around  industrial  plants,  and  atmospheric  trace  chemicals 
in  the  stratosphere.  Lidar  has  also  been  used  to  measure  the  velocity 
and  direction  of  winds  near  storms  and  airports  and  to  track  the 
global  circulation  of  volcanic  ash  emitted  into  the  atmosphere  after 
recent  eruptions,  such  as  those  at  Mount  St.  Helens  and  El  Chichon. 
Under  optimal  conditions  lidar  can  be  extremely  sensitive  An 
example  is  the  ground-based  laser  remote  sensing  of  sodium  and 
lithium  atoms  in  the  stratosphere  at  ranges  greater  than  90  km  and 
in  concentrations  as  low  as  a  few  atoms  per  cubic  centimeter  (5). 
More  commonly,  detection  ranges  arc  on  the  order  of  a  few  hundred 
meters  to  several  kilometers  and  concentration  levels  on  the  order  of 
parts  per  million  to  parts  per  billion 
The  use  of  optical  backscattcr  to  measure  properties  of  the 
atmosphere  is  not  new,  extensive  experiments  having  been  conduct¬ 
ed  in  the  early  1900’s  with  large  searchlights.  The  field  of  optical 
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remote  sensing  was  greatly  advanced  by  the  laser,  which  offers 
several  improvements  over  conventional  light  sources.  These  include 
narrow  spectral  width  (<0.01  nm),  a  frequency  or  color  that  is  often 
tunable,  and  high  peak  power  (>106  W)  available  in  a  short  pulse 
(<1  jxsec)  and  in  a  narrow  beam  (<10  cm  in  diameter).  These 
attributes  make  the  laser  an  ideal  spectroscopic  probe  of  the 
atmosphere.  In  this  regard,  a  lidar  system  may  be  thought  of  as  an 
“active”  remote  sensing  system  since  it  can  illuminate  the  target 
region,  in  contrast  to  a  “passive”  optical  sensor  which  deteas 
ambient  light  or  thermal  emission  from  the  target. 

The  applications  of  lidar  systems  for  the  remote  sensing  of 
atmospheric  properties  were  appreciated  soon  after  the  discovery  of 
lasers  in  the  early  1960’s.  Early  lidar  measurements  were  made  in 
1962  by  Fiocco  and  Smullin  (7)  who  bounced  a  laser  beam  off  the 
moon  and  who  also  investigated  the  turbid  layers  in  the  upper 
atmosphere;  in  1963  by  Ligda  (S)  who  used  a  rubv  laser  to  obtain 
the  first  lidar  measurements  of  cloud  heights  and  tropospheric 
aerosols,  and  in  1964  by  Schotiand  (9)  who  used  a  temperature- 
tuned  ruby  laser  to  detea  water  vapor  in  the  atmosphere.  Progress 
has  been  continuous  since  that  time,  but  the  discovers'  of  different 
laser  sources  in  the  past  decade,  coupled  to  improvements  in  optical 
instrumentation  and  data  processing,  has  been  responsible  for  the 
recent  surge  in  the  number  of  laser  remote  sensing  systems.  This 
improved  capability  has  been  accompanied  by  an  increased  aware¬ 
ness  of  the  need  to  monitor  the  impact  of  natural  and  anthropogenic 
influences  on  the  environment. 

Remote  sensing  of  the  atmosphere  by  optical  techniques  can  be 
accomplished  in  several  ways.  One  technique  involves  measurement 
of  the  absorption  spearum  of  the  atmosphere  over  a  long  path 
separating  a  spcaroscopic  optical  source  and  a  detector  (5).  Another 
long-path  absorption  technique  uses  a  configuration,  described  by 
HinkJev  and  Kelley  (10),  in  which  a  tunable  laser  source  and 
dctcaor  arc  located  together  and  a  rctroreflccting  mirror  is  placed  at 
a  distance  of  several  hundred  meters;  such  a  system  is  useful  when 
the  laser  source  is  weak,  since  the  retroreflector  gready  enhances  the 
returned  radiation.  In  this  article  we  will  primarily  be  concerned 
with  a  pulsed  lidar  system  where  the  laser  and  detector  arc  located 
together  and  no  retroreflector  is  used  as  the  target;  in  this  case  the 
returned  laser  radiation  is  due  to  backscatter  from  aerosols  or  dust  in 
the  atmosphere  or  a  topographic  target  such  as  a  hill  or  trees. 
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Optical  Requirements  for  Laser  Remote 
Sensing 

The  physics  and  optical  technology  associated  with  atmospheric 
laser  remote  sensing  encompass  atmospheric  optical  propagation, 
atomic  and  molecular  spectroscopy,  digital  signal  processing,  and 
laser  technology.  In  this  section  we  discuss  the  optical  requirements 
for  iidar  and,  in  particular,  review  (i)  the  presence  of  transmission 
“windows”  in  the  atmosphere  which  permit  the  laser  beam  to 
propagate  through  the  atmosphere  with  little  attenuation,  (ii)  the 
identification  and  relative  strength  of  the  various  optical  interactions 
between  an  atmospheric  species  and  a  laser  beam,  and  (iii)  the 
properties  of  suitable  lasers. 

Atmospheric  transmission  windows.  For  a  laser  beam  to  detea  a 
species  at  a  long  distance,  the  beam  must  not  be  appreciably 
attenuated  by  the  intervening  atmosphere.  The  output  wavelength 
of  the  laser  must  therefore  lie  in  a  spectral  transmission  “window”  of 
the  atmosphere.  These  windows  ?re  shown  as  clear  areas  in  Fig.  2, 
which  shows  a  low-resolution  transmission  spectrum  of  the  atmo¬ 
sphere  over  the  spectral  region  from  the  ultraviolet  to  the  mid- 
infrared  (15  pm)  for  a  path  of  300  m  .(11).  The  atmosphere  has 
regions  of  strong  absorption  as  well  as  transmission  windows;  the 
absorption  regions  are  due  primarily  to  oxygen,  carbon  dioxide,  and 
water  vapor,  and  exhibit  finer  detail  than  that  shown  in  the  figure. 
The  most  useful  transparent  spectral  ranges  of  the  atmosphere  are 
the  visible  (0.4  to  0.7  pm),  near  infrared  (0  7  to  1.5  pm),  and  3-  to 
5-(xm  and  9-  to  13-pm  regions.  Within  these  spectral  regions,  the 
laser  is  not  appreciably  attenuated  except  by  clouds,  and  thus  remote 
sensing  over  long  distances  may  be  achieved. 

Optical  interaction  of  lasers  with  atmospheric  species.  The  primary 
influence  of  the  atmosphere  on  a  low- power  laser  beam  is  through 
scattering  and  absorption.  Both  processes  cause  an  attenuation  of 
the  beam  according  to  Bougueris  or  Beer’s  law 

i  =  (1) 

where  7  is  the  intensity  of  the  optical  beam  after  transmission  over  a 
distance  it,  a  is  the  atmospheric  extinction  coefficient,  and  70  is  the 
initial  intensity  of  the  beam.  It  is  possible  to  express  a  as  a  sum  of 
terms 


O  —  OlRay  +  OtMie  aRam*n  +  ®ab$ 


(2) 


where  aRay,  aM,e,  and  aRair)an  arc  the  extinaion  coefficients  related 
to  Rayleigh,  Mic,  and  Raman  scattering,  respectively,  and  aabs  is  the 
molecular  absorption  coefficient  (3).  Rayleigh  scattering  is  due  to 
parades  in  the  atmosphere,  such  as  molecules  or  fine  dust,  that  are 
much  smaller  than  the  optical  wavelength,  X.,  of  the  laser.  Mie 
scattering  is  associated  with  larger  parades  such  as  aerosols  whose 
size  is  on  the  order  of  X.  Rayleigh  and  Mic  processes  arc  elastic 
scattering:  the  scattered  light  is  the  same  wavelength  (color)  as  the 
incident  laser  beam  Raman  scattering  is  an  inelastic  interaction  of 
the  optical  beam  involving  excitation  of  the  energy  levels  of  a 
molecule  and  rcradiation  at  a  different  wavelength  Absorption  of 
the  laser  beam  is  a  resonant  interaction  (direa  absorption)  leading 
to  excitation  of  the  molecule,  followed  possibly  by  fluorescence 
The  sensitivity  of  laser  remote  sensing  is  related  to  the  relative 
strength  of  these  optical  processes  through  the  “lidar  equation.”  The 
Iidar  equation  (2,  3)  is  used  to  predia  the  power  of  the  backscat- 
tered  signal,  Pr,  and  may  be  given  in  the  form 


Pr  = 


PipAKe-2aJI 

ir*2 


(3) 


where  P,  is  the  transmitted  laser  power,  R  is  the  range  to  the  target, 
p  is  the  effective  reflectivity  of  the  target  (assumed  to  be  larger  than 
the  laser  beam),  which  backscattcrs  the  radiation  into  v  stcradians, 


Detector 


Fig.  i  .  Schematic  of  a  Iidar  used  for  laser  remote  sensing  of  the  atmosphere. 
The  backscartercd  laser  radiation  may  contain  spectral  information  on  the 
composition  and  physical  condition  of  the  cloud  or  the  intervening  atmo¬ 
sphere. 


A  is  the  area  of  the  telescope,  K  is  the  optical  efficiency  of  the  system, 
and  cxp(-2aft)  is  the  two-way  atmospheric  attenuation  of  the 
optical  beam.  The  more  general  form  of  Eq  3  contains  integrals 
over  range,  time,  and  spectral  wavelength  (<5).  The  target  can  be 
either  a  solid  refleaor,  such  as  a  building  or  trees,  or  a  “distributed” 
refleaor  consisting  of  molecules,  dust,  or  clouds  in  the  atmosphere. 
In  the  latter  case,  the  target  reflectivity  term,  p,  is  due  to  Mie, 
Rayleigh,  Raman,  or  fluorescence  processes  and  is  thus  also  related 
to  the  associated  a  values. 

The  approximate  strength  of  these  processes  is  given  in  Table  1, 
which  lists  typical  maximum  values  of  p  and  a  for  some  of  the  most 
important  optical  interactions  (3,  6)  The  contribution  of  the 
scattering  processes  to  both  a  and  p  arises  because  scattering  reduces 
the  intensity  of  the  laser  beam,  contributing  to  extinction,  whereas 
the  portion  of  the  beam  that  is  scattered  by  180°  contributes  to 
the  deteaed  signal,  thereby  serving  effcaivcly  as  a  refleaor.  Absorp¬ 
tion  contributes  onlv  to  the  exnnaion  coefficient,  with  no  intrinsic 
contribution  to  p.  In  Table  1,  values  arc  given  for  p  in  conjunction 
with  absorption,  and  refer  to  effective  reflectivities  arising  from 


Table  1.  Optical  interactions  for  laser  remote  sensing  of  the  atmosphere; 
N.A.,  not  applicable. 


Optical 

interaction 

Spectral 

range 

(M-m) 

Scartcrcr 

Attenua¬ 
tion*  a 
(km’1) 

Effcctivef 
reflectivity  p 

Absorption 

UV 

Atoms, 

107 

10~2 

(resonant) 

I 

molecules 

106 

10"3  to  10“l 

10 

103 

10‘7  to  10'4 

Scattering  (non¬ 

resonant) 

Rrylcigh 

UV 

o2,  n2 

10'1 

lO-3 

1 

10'3 

10  3 

10 

10"7 

10"9 

Mic 

UV 

Aerosols, 

10'1 

nr3 

1 

fine  dust 

10~' 

10’4 

10 

10-2 

10“3 

Raman 

UV 

Molecules 

10  5 

10~7 

1 

10-7 

10-9 

10 

10-'1 

10-'5 

B^ckscatter 

UV-1R 

Buildings, 

N.A. 

10  2  to  1 

by  hard 

trees, 

targets 

dense  clouds 

♦Attenuation  calculated  for  concentration  of  1  atm.  +  Reflectivity  calculated  for 
range  resolution  of  100  m,  1  atm  concentration,  and  isotropic  scattering. 
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fluorescence  emission  of  the  atom  or  molecule  after  it  has  ab¬ 
sorbed  the  incident  light  (12).  Conversely,  a  topographic  target  docs 
not  contribute  to  signal  extinction;  instead,  the  180°  scattered 
portion  of  the  incident  beam  corresponds  to  the  effective  target 
reflectivity. 

In  general,  the  value  of  p  or  a  may  be  a  complicated  function  of 
’laser  lincwidth,  wavelength,  or  spectral  filtering  in  the  detection 
process,  and  may  exhibit  strong  resonances.  The  values  presented  in 
Table  1  are  appropriate  to  the  strong  resonance  case  and  are 
appropriate  only  if  the  laser  spectral  requirements  arc  met.  The  large 
values  for  a  shown  for  absorption  are  for  the  case  where  the  laser 
wavelength  is  coincident  with  an  absorption  line  of  the  molecule; 
off- resonant  values  of  a  are  orders  of  magnitude  smaller.  The 
backscatter  and  absorption  coefficients,  in  general,  depend  on 
species  concentration,  temperature,  pressure,  and  possibly  velocity, 
so  that  lidar  can  measure  one  or  more  of  these  parameters.  Measure¬ 
ments  for  ail  lidar  systems  are  based  on  changes  in  the  returned 
signal,  Pn  as  a  function  of  wavelength  or  range,  and  thus  deduce 
information  about  the  physical  state  of  the  atmosphere  through  the 
effect  of  a  or  p. 

A  wide  range  of  values  occur  for  ot  and  p  (Table  1).  Since  all 
optical  processes  shown  can  occur  simultaneously,  it  is  important  to 
design  the  lidar  system,  through  judicious  choice  of  laser  wave¬ 
length,  spectral  filtering,  or  geometry,  so  that  the  measurement  is 
most  sensitive  to  only  one  optical  interaction  In  some  lidar  designs, 
a  and  p  depend  on  different  optical  interactions;  in  other  words,  the 
attenuation  and  backscatter  arise  from  different  processes. 

Lidar  systems  that  measure  a  particular  optical  interaction  have 
been  classified  as  follows:  (i)  atmospheric  backscatter  lidar,  (ii) 
differential -absorption  lidar  (DIAL),  (iii)  fluorescence  lidar,  (iv) 
Raman  lidar,  and  (v)  Doppler  lidar.  These  classifications  refer  to  the 
optical  interaction  to  which  the  system  has  been  designed  to  be  most 
sensitive.  Atmospheric  backscatter  lidar  transmits  one  laser  wave¬ 
length  and  detects  changes  in  the  backscatter  due  to  the  aerosols  or 
dust  in  the  atmosphere.  DIAL  measures  the  concentration  of  a 
molecular  species  in  the  atmosphere  by  transmitting  twro  wave¬ 
lengths,  only  one  of  which  is  absorbed,  and  detecting  the  difference 
in  the  intensity  of  the  returns  at  the  two  wavelengths;  the  backscat- 
tcr  in  DIAL  may  be  from  a  hard  target,  or  aerosol  and  dust  returns 
as  in  atmospheric  backscatter  lidar  Fluorescence  lidar  uses  two 
wavelengths  as  in  a  DIAL  system  and  in  addition  uses  spectromctric 
techniques  to  separate  the  wavelength -shifted  fluorescence  signal 
from  the  strong  Rayleigh  backscatter  in  the  atmosphere.  Raman  and 
Doppler  lidar  uses  a  single -wave length  laser  and  sophisticated 
spectromctric  detection  techniques  to  spectrally  resolve  the  wave¬ 
length-shifted  signal  from  the  strong  background  due  to  Rayleigh  or 
Mic  scattering.  A  Raman  measurement  often  detects  a  wide  spectral 
range  (~I0  nm)  containing  a  number  of  vibrational -rotational 
Raman-shifted  lines,  whereas  a  Doppler  measurement  detects  only  a 
very  narrow  spectral  range  (~10~5  nm)  that  encompasses  the 
Doppler-shifted  backscatter  lidar  return.  Because  DIAL,  fluores¬ 
cence,  and  Raman  techniques  use  several  wavelengths  that  arc  cither 
transmitted  by  the  laser  or  detected  by  the  lidar  system,  they  provide 
measurements  of  the  absolute  concentration  of  the  detected  species. 
In  contrast,  atmospheric  backscatter  lidar  uses  a  single-wavelength 
laser  and  single-wavelength  detection  and  is  best  suited  for  lidar 
measurements  that  determine  relative  changes  in  the  density  of  a 
known  atmospheric  species;  these  relative  density  measurements 
may,  however,  be  referenced  to  an  absolute  value  if  the  lidar  is 
calibrated  with  the  use  of  a  standard  target. 

Laser  sources.  The  utility  of  laser  remote  sensing  is  dependent  on 
several  factors,  including  the  availability  of  lasers  with  sufficient 
output  power  and  wavelength  characteristics  to  provide  an  adequate 
lidar  signal.  Although  many  lasers  provide  high  power,  they  often 
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Fig.  2.  Optical  transmission  of  the  atmosphere  depicting  the  presence  of 
strong  absorption  bands  and  transmission  “windows”  in  die  atmosphere. 
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Fig.  3.  Typical  values  of  laser  pulse  energy  as  a  funcuon  of  laser  wavelength 
tuning  range  showing  the  gaps  in  wavelength  coverage  of  currently  available 
lasers;  (  +  )  line  tunable;  (*)  continuously  tunable. 


lack  spectral  coverage  or  tunability,  preventing  a  determination  of 
the  resonant  v  alues  of  a  or  p  (Table  1)  for  the  particular  atmospheric 
species  or  for  the  optical  parameter  of  interest. 

This  may  be  better  understood  on  the  basis  of  Fig.  3,  which 
presents  the  output  energy  per  pulse  now  available  from  commer¬ 
cially  available  pulsed  lasers  of  moderate  size  as  a  function  of  their 
wavelength  coverage  (13,  14).  The  wavelength  tuning  range  of 
high-power  lasers  is  not  complete  across  the  spectral  range  shown 
and  is  particularly  lacking  in  the  “eye-safe”  spectral  region  (X  >  1.4 
pm)  For  remote  sensing  in  the  infrared  (IR),  a  minimum  transmit¬ 
ted  laser  pulse  energy  on  the  order  of  1  to  10  mj  is  required  for  the 
detection  of  lidar  signals  from  a  hard  target  at  a  range  of  a  few 
kilometers,  and  about  1  J  for  range- resolved  returns  from  atmo¬ 
spheric  aerosols.  These  values  may  be  be  lowered  by  one  to  two 
orders  of  magnitude  if  heterodyne  detection  can  be  used  to  detect 
the  IR  lidar  return  or  if  the  remote  sensing  can  be  conducted  in  the 
visible  or  ultraviolet  (UV).  The  detection  of  lidar  atmospheric 
backscatter  is  much  more  sensitive  in  the  latter  spectra!  regions 
because  detectors  less  susceptible  to  noise  (such  as  photomultiplier 
tubes)  arc  available  and  the  atmospheric  backscatter  is  improved  in 
the  UV-visible  as  compared  to  the  IR. 

The  characteristics  shown  in  Fig.  3  demonstrate  why  most  laser 
remote  sensing  experiments  have  used  ruby,  dye,  neodymium: 
yttrium-aluminum-gamet  (Nd.YAG),  and  C02  lasers.  These  lasers 
offer  sufficient  energy  for  the  detection  of  atmospheric  species  at 
moderate  ranges,  operate  at  wavelengths  that  arc  not  appreciably 
attenuated  by  the  atmosphere,  and  occur  in  a  spectral  range  where 
several  atmospheric  constituents  have  reasonably  strong  optical 
interactions.  Some  of  the  lasers  shown,  such  as  Nd.YAG,  CO,  and 
C02,  arc  line  tunable  in  that  they  operate  only  on  a  limited  number 
of  discrete  lines  or  wavelengths.  The  other  lasers  are  continuously 
tunable  within  the  spectral  range  shown  in  Fig.  3.  Single-wave¬ 
length  lasers  are  suitable  for  atmospheric  backscatter,  Raman,  and 
Doppler  lidar  measurements  since  the  spectral  information  is  ob¬ 
tained  during  the  detection  process  with  a  spectrometer  or  spectral 
filter.  Lasers  for  DIAL  or  fluorescence  lidar  must  be  tunable, 
however,  so  that  the  laser  wavelength  can  match  the  resonance 
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Fig.  4.  Airborne  lidar  returns  as  a  function  of  altitude  and  longitude  near  the  coast  of  South  Carolina  on  the  evening  of  27  May  1980  show  ing  the  in¬ 
creased  lidar  backscacter  in  a  layer  near  14  to  16  km  due  to  volcanic  material  injected  into  the  atmosphere  by  Mount  St.  Helens  volcano.  [Adapted  from 
(16)} 
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Fig.  5.  Airborne  lidar  data  near  Los  Angeles  showing  the  vertical  density  profile  of  the  boundary  layer  and  continental  flow  of  urban  haz*.  [Adapted  from 


(17)} 

wavelength  of  the  optical  absorption  lines  of  the  molecules  or  atoms 
being  detected.  This  is  most  important  for  the  remote  sensing  of 
molecules  or  pollutants  in  the  atmosphere  because  most  of  the 
absorption  lines  of  many  important  organic  species  occur  in  the 
middle  IR  (2  to  5  fxm)  and  this  “fingerprint”  spectral  region  does 
not  coincide  with  many  high -power,  tunable  laser  sources.  Because 
of  these  considerations,  one  of  the  most  important  goals  in  the 
future  application  of  laser  remote  sensing  is  the  development  of 
new,  high-power  tunable  lasers  to  cover  the  spectral  gaps  shown  in 

Fig.  3- 


Lidar  Atmospheric  Measurements 

Each  of  the  five  lidar  types  can  be  used  for  the  remote  sensing  of  a 
selected  variety  of  atmospheric  constituents  or  parameters.  Table  2 
lists  some  of  the  constituents  and  parameters  that  have  been 
measured  with  lidar.  More  detail  is  available  in  the  technical  digests 
of  several  recent  lidar  conferences  (15). 

Atmospheric  bacbscatter  lidar  Atmospheric  backscattcr  lidar  is  the 
most  common  type  of  lidar  system  and  consists  of  a  nontunablc, 
high- power,  pulsed  laser  that  is  used  to  measure  changes  in  the 
density  or  backscattcr  coefficient  of  the  atmosphere  that  arc  due  to 
particulates,  such  as  aerosols,  clouds,  dust,  or  plumes.  These  atmo¬ 
spheric  constituents  have  relatively  large  optical  scattering  cross 
sections  and  arc  relatively  easy  to  detect  (Table  1).  Such  lidars  have 
been  useful  in  tracking  the  turbid  effluent  and  gas  plumes  from 
industrial  smokestacks  and  have  also  been  used  to  map  out  rain, 
sr.ow,  ice  crystals,  and  dense  clouds  in  the  atmosphere  (1).  Precision 
measurements  of  the  slight  changes  in  the  vertical  density  of  the 
atmosphere  have  also  been  made  with  lidar,  from  which  one  can 
determine  the  pressure  and  temperature  of  the  atmosphere  as  a 
function  of  altitude  (15). 

Recently,  McCormick  used  a  lidar  system  to  track  and  map  out 
the  density  of  volcanic  ash  emitted  into  the  atmosphere  by  Mount 
St.  Helens  and  El  Chichon  (16).  Figure  4  presents  data  obtained  by 
McCormick  using  an  airborne  lidar  system  incorporating  a  high- 
power,  frequency-doubled  Nd:YAG  laser  operating  near  0.53  pm 
(16).  A  thick  layer  of  volcanic  ash  was  detected  at  an  altitude  near  16 


km  and  was  found  to  vary  in  density  and  stratification  as  a  function 
of  distance  By  tracking  the  cloud  of  volcanic  ash,  McCormick  was 
able  to  use  the  ash  as  a  tracer  for  global  wind  and  circulation 
patterns.  Such  data,  coupled  with  IR  satellite  data,  showed  the 
kinetics  of  global  atmospheric  circulation  and  tracked  the  slow 
global  mixing  of  the  volcanic  ash  between  the  Northern  Hemisphere 
and  the  Southern  Hemisphere  that  occurred  the  following  year. 

Atmospheric  aerosols  have  also  been  studied  by  lidar  in  both  the 
troposphere  and  the  stratosphere.  Lidar  systems  with  beam-steering 
capabilities  have  been  used  by  Uthe  et  at.  for  three-dimensional 
mapping  of  the  flow  of  dense  tropospheric  aerosols  (haze)  near  the 
ocean-land  boundary  and  urban  centers  (17).  An  example  of  such 
measurements  obtained  near  Los  Angeles,  showing  the  flow  pat¬ 
terns  of  urban  aerosols,  is  given  in  Fig.  5  (17)  Such  data  have 
greativ  increased  our  knowledge  of  the  interplay  of  air  circulation 
and  pollution  emission  and  arc  being  used  to  predict  the  impact  of 
emissions  near  urban  centers. 

Atmospheric  backscattcr  lidar  has  been  instrumental  in  investiga¬ 
tions  of  the  propagation  of  a  laser  beam  through  the  atmosphere. 
These  investigations  have  included  the  study  of  complicated  scatter¬ 
ing  processes  (1)  present  under  normal  atmospheric  conditions  due 
to  the  statistical  distribution  and  variability  of' aerosols  and  scatter¬ 
ing  particles  (18).  In  addition,  the  effect  of  multiple  scattering  has 
been  investigated  by  lidar;  multiple  scattering  is  due  to  optically 
thick  clouds  that  can  drastically  alter  the  transmission  properties  of  a 
laser  beam  from  that  described  by  Eq.  I  (19). 

Differential-absorption  lidar.  Most  lidar  measurements  that  deter¬ 
mine  the  concentration  of  gaseous  species  in  the  atmosphere  are 
presently  made  with  the  DIAL  technique  (20).  This  technique  uses 
two  or  more  wavelengths  to  measure  the  difference  in  the  absorp¬ 
tion  of  the  lidar  signal  as  the  laser  frequency  is  varied  between  a 
wavelength  that  is  absorbed  by  molecules  in  the  atmosphere  and  a 
wavelength  that  is  not  absorbed.  The  method  is  sensitive  because  it 
takes  advantage  of  the  relatively  large  value  of  the  “on-resonance” 
absorption  coefficient  a  (Table  I). 

In  principle,  the  concentration  of  an  atmospheric  species  may  be 
obtained  directly  from  Eqs.  2  and  3  on  the  basis  of  measurements  at 
a  single  wavelength.  However,  a  concentration  determination  based 
on  Eq  3  and  a  single-wavelength  measurement  requires  an  accurate 
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knowledge  of  K,  p,  a,  and  the  geometric  factors  R  and  A.  In 
general,  these  parameters  arc  either  unknown  or  known  with  poor 
accuracy 

The  DIAL  technique  attempts  to  overcome  this  problem  by 
measuring  the  Ildar  signals  at  two  wavelengths,  \  and  X\  The  laser 
wavelength  \  is  chosen  to  coincide  with  a  wavelength  strongly 
absorbed  by  the  molecule  being  investigated,  while  at  X'  there  is 
little  or  no  absorption.  The  use  of  two  wavelengths  results  in  two 
equations  of  the  form  of  Eq.  3,  which  arc  then  solved  simultaneous¬ 
ly  for  Na,  the  path- averaged  concentration  of  the  absorbing  mole¬ 
cule  to  yield 


1 


2(<ra'  -  aa)it 


KfMf) 


+  2(n  -  a')R 


(4) 


where  cra  is  the  absorption  cross  section  of  the  absorbing  molecule, 
P  is  the  backscattcrcd  signal  power  normalized  to  the  transmitted 
power  ( P  =  iyTt),  and  the  primed  and  unprimed  parameters  refer 
to  values  at  the  low-  and  high-absorption  wavelengths  X'  and  X, 
respectively  (20).  In  deriving  Eq.  4  from  Eq,  3,  the  e~2aJt  term  in 
Eq.  3  has  been  replaced  by  where  a  now  represents  the 

assumed  known  total  attenuation  of  the  atmosphere  except  that  due 
to  absorption,  and  Na<ra  is  the  attenuation  due  to  absorption  by  the 
molecule  being  investigated. 

The  advantage  of  this  rwo-ffcqucncy  (DIAL)  approach  is  that 
only  differences  in  the  various  parameters  need  be  considered.  For 
this  reason  the  DIAL  technique  is  more  commonly  used  for  species 
concentration  determination  than  single-frequency  absorption  mea¬ 
surements,  despite  the  added  complexity  of  a  dual- frequency'  lidar 
system.  Further  simplification  can  be  achieved  when  the  use  of 
closely  spaced  frequencies  in  a  given  measurement  results  in  differ¬ 
ences  of  target  reflectivity  (Ap)  and  atmospheric  background  extinc¬ 
tion  (Aa)  that  are  negligibly  small.  In  that  case,  Eq.  4  reduces  to 


ln(iyPa') 
2(aa  -  cra')il 


(5) 


The  form  given  in  Eq  5,  which  greatly  simplifies  the  relation 
between  the  lidar  return  ratio  and  species  concentration,  is  frequent¬ 
ly  assumed  in  DIAL  measurements,  but  Eq.  4  must  be  used  in  eases 
where  either  Ap  or  Aa  is  nonnegligiblc.  Equation  5  is  valid  for 
column-content  or  path-averaged  DIAL  measurements  and  is  slight¬ 
ly  changed  for  range- resolved  measurements. 

A  large  number  of  molecular  species  have  been  studied  with  the 
DIAL  technique,  and  the  results  of  many  of  these  experiments  arc 
given  in  the  proceedings  of  recent  conferences  (75),  including 
extensive  review  articles  on  DIAL  by  T.  Wilkcrson  and  by  E. 
Browcll.  The  lidar  systems  used  in  these  DIAL  measurements  span 
the  frequency  range  from  the  UV  to  the  IR  and  cover  both  rangc- 
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Fig.  6.  C02  DIAL  laser  remote  sensing  of  atmospheric  CO  for  a  path  over  a 
traffic  roadway.  Roadway  range,  480  m;  lidar  range,  500  m. 

resolved  (RR)  and  path-averaged  or  column-content  (CC)  measure¬ 
ments. 

Most  DIAL  experiments  to  date  have  used  lidar  signals  backscat¬ 
tcrcd  from  cither  a  retroreflcctor  or  topographic  targets  (trees, 
buildings)  to  measure  the  column  content  of  the  atmospheric  species 
between  the  laser  transmitter  and  the  target.  Furthermore,  mosr  of 
these  experiments  have  been  earned  out  in  the  IR,  where  almost  all 
molecules  of  interest  have  extensive  absorption  bands.  The  mole¬ 
cules  investigated  in  CC  DIAL  experiments  include  SO2,  NH3,  O3, 
CO,  C02,  HCI,  H20,  NO,  N2H4,  N20,  and  SF6  (4-6,  21).  The 
sensitivity  of  many  of  these  CC  DIAL  measurements  is  on  the  order 
of  a  few  parts  per  billion  to  a  few  parts  per  million,  with  detection 
ranges  on  the  order  of  Several  kilometers.  As  an  example,  Fig.  6 
shows  time-averaged  CC  DIAL  data  obtained  for  the  remote 
sensing  of  CO  emitted  by  automobiles  at  a  range  of  500  m  with  a 
tunable  C02  laser  DIAL  system  (22)  .  The  increase  in  CO  concentra¬ 
tion  dunng  the  morning,  lunch,  and  evening  rush  hours  is  evident. 
The  sensitivity  of  the  DIAL  measurement  was  approximately  10 
ppb  and  the  temporal  resolution  was  approximately  1  second;  the 
CO  emission  from  individual  vehicles  was  easily  observed. 

Range -resolved  DIAL  measurements  provide  a  continuous  signal 
as  a  function  of  range,  based  on  the  backscarter  from  atmospheric 
aerosols  and  particulates.  As  such,  RR  DIAL  measurements  have  a 
significant  advantage  over  CC  DIAL  measurements  in  that  they  can 
provide  a  three-dimensional  map  of  the  measured  species.  The  major 
difficulty  with  RR  DIAL  is  the  relatively  small  backscarter  coeffi¬ 
cient  for  atmospheric  aerosols  and  particulates  compared  to  that  for 


Table  2.  Selected  list  of  atmospheric  constituents  and  parameters  measured  by  lidar.  Accuracy  and  ranges  given  arc  typical  values  uid  depend  upon  individual 
lidar  measurements. 


Constituent 
or  parameter 

Laser  type 

Accuracy 

Lidar  type 

Range 

(km) 

Dust,  clouds,  volcanic  ash,  smoke  plumes 

Ruby,  Nd:YAG 

I  to  10% 

Atmospheric  backscarter 

10  to  50 

h2o,  0,,  so2,  no,  no2,  n2o,  c2h<, 

CH4,  HCI,  CO,  Hg,  SF6,  NH3,  C02 

Dye.  C02,  OPO,* 
cxcimcr,  Co:MgF2 

Variable,  I  ppb  to 

100  ppm 

DIAL,  Raman 

1  to  5 

OH,  Na,  K,  Li,  Ca,  Ca+ 

Dye 

I02  to  107  atoms  per 
cubic  centimeter 

Fluorescence 

1  to  90 

Temperature,  pressure 

Dye,  Nd  YAG 

1  K,  5  mbar 

DIAL,  Raman 

1  to  30 

Wind  speed 

co2 

0.5  m/scc 

Doppler 

15 

*OPO,  optical  parametric  oscillator. 
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Fig.  7.  Illustrations  of  the  laser  remote  sensing  of  the  concentration  of 
ethylene  near  the  distillation  towers  of  a  refinery  (Left)  Isobars  of  concen¬ 
tration  in  increments  of  20  ppb.  (Right)  Corresponding  three-dimensional 
plot  showing  a  peak  concentration  of  approximately  150  ppb.  [Adapted 
from  (27)] 

a  hard  target.  Because  of  the  small  backscarter  coefficients,  most  RR 
DIAL  measurements  have  used  relatively  high-power  lasers,  or 
sophisticated  optical  detection  techniques  such  as  photon  counting 
or  heterodyne  detection  to  increase  sensitivity  to  the  backscartcrcd 
signal.  In  spite  of  these  difficulties,  several  RR  DIAL  measurements 
of  atmospheric  molecules  have  been  made.  These  measurements 
include  the  detection  and  mapping  of  SO2  and  NO2  in  a  power  plant 
plume  (5),  SO2  and  O3  near  urban  centers  (23),  and  N02  emitted  by 
a  chemical  factory  (24).  Other  examples  arc  atmospheric  tempera¬ 
ture  and  pressure  measurements  based  on  absorption  lines  of  oxygen 
near  760  run  (25)  and  precision  measurements  of  O3  in  the  upper 
atmosphere  (26)  These  O3  measurements  arc  particularly  important 
in  light  of  the  recent  discovery  of  a  “hole”  in  the  atmospheric  O3 
layer  near  Antarctica  and  the  anticipated  use  of  a  similar  DIAL 
system  bv  NASA  to  study  this  phenomenon 

Finallv,  as  an  example  of  RR  DIAL  measurements,  Fig.  7  shows 
data  obtained  by  Rothe  using  a  CO2  laser  DIAL  system  (27)  which 
measured  the  concentration  of  ethylene  (C2H4)  over  the  area  of  a 
refinery.  This  figure  illustrates  the  three-dimensional  mapping  capa¬ 
bility  of  such  a  system. 

Fluorescence  lidar.  In  a  fluorescence  lidar  system,  the  laser  is  tuned 
to  an  absorption  line  of  the  species  that  is  to  be  measured  and  the 
reradiated  fluorescence  is  detected  by  selective  spectral  filtering  of 
the  returned  radiation.  The  fluorescence  radiation  may  be  at  the 
same  wavelength  as  the  excitation  wavelength  or  may  be  red-shifted 
(have  a  longer  wavelength).  In  the  lartcr  ease,  a  spectrometer  or 
narrow-band  interference  filter  is  usually  used  to  reject  the  strong 
Rayleigh  backscarter  and  thus  detea  only  the  fluorescence  signal. 
The  backscarter  coefficient  for  fluorescence  is  much  greater  in  the 
UV  than  in  the  IR  (Table  1);  this  is  due  to  the  combined  effects  of 
absorption  cross  seaion  (greater  in  the  UV  than  the  IR)  and 
radiative  lifetime  (longer  in  the  I R  than  in  the  UV,  which  permits 
increased  de-excitation  due  to  collisions  with  ambient  molecules). 

Several  faaors  that  limit  the  application  of  fluorescence  lidar  for 
remote  sensing  include  problems  involving  detcaor  sensitivity 
coupled  with  solar  background  radiation,  which  generally  reduces 
fluorescence  measurements  for  remote  sensing  to  nighttime  investi¬ 
gations  and  to  wavelengths  shorter  than  1  p,m,  where  photomulti¬ 
plier  dctcaion  may  be  used.  In  addition,  at  tropospheric  pressures 
an  excited  molecule  will  normally  undergo  several  collisions  with 
other  molecules  before  fluorescent  emission.  This  quenching  of  the 
excited  molecule  severely  reduces  the  fluorescent  output,  rendering 
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it  extremely  ineffiaent  for  most  remote  sensing  applications  in  the 
troposphere  or  even  in  the  lower  stratosphere. 

In  spite  of  these  difficulties,  several  successful  measurements  of 
trace  atmospheric  species  have  been  made  with  this  technique 
Ground-based  visible  fluorescence  lidar  systems  have  been  devel¬ 
oped  by  Megie  (28)  and  by  Gardner  and  his  colleagues  (29)  and 
have  been  used  to  study  alkali  metal  (sodium,  potassium,  lithium, 
calcium,  Ca+)  concentration  profiles  at  altitudes  of  80  to  100  km  (5, 
28 t  29).  Peak  sodium  concentrations  of  about  I03  to  104  atoms  per 
cubic  centimeter  have  been  measured  in  the  upper  atmosphere, 
potassium  concentrations,  about  an  order  of  magnitude  lower,  and 
lithium  atoms,  with  peak  concentrations  of  approximately  I  atom 
per  cubic  centimeter  at  altitudes  near  100  km,  have  been  observed. 
The  return  signal  from  sodium  atoms  is  sufficient  to  serve  as  a  tracer 
to  study  temperature,  eddy  diffusion,  and  the  dynamics  of  the 
wavelike  features  believed  to  be  due  to  internal  gravity  waves  and 
tidal  effects  in  the  upper  atmosphere. 

Another  significant  application  of  fluorescence  lidar  involves  the 
study  of  the  hydroxyl  free  radical,  OH.  Although  OH  occurs  as  a 
trace  species  at  very  low  concentrations  (on  the  order  of  0. 1  part  per 
trillion),  it  plays  an  important  catalytic  role  in  various  atmospheric 
chemistry  processes  and,  together  with  chlorine  and  nitrogen  ox¬ 
ides,  is  involved  m  the  O3  destruction  cycle  in  the  atmosphere. 
Airborne  measurements  of  OH  have  been  made  with  a  laser 
operating  at  282  nm  (30).  Concentrations  on  the  order  of  IO6  to  107 
OH  per  cubic  centimeter  were  observed,  corresponding  to  approxi¬ 
mately  one  part  in  1013.  These  measurements  complement  the 
extensive  investigation  of  OH  in  the  stratosphere  in  which  a  laser 
induced  fluorescence  spectrometer  was  mounted  on  a  balloon-borne 
platform  (31). 

Raman  scattering  lidar.  The  use  of  Raman  scattering  methods  for 
remote  sensing  is  severely  limited  by  the  small  optical  interaction 
strength  for  Raman  scattering.  As  a  result,  this  technique  is  generally 
used  with  high-energy  pulsed  lasers,  is  rcstnacd  to  thv  UV  or  visible 
regions  of  the  spectrum  so  that  sensitive  photomultiplier  tubes  can 
be  used  for  detection,  and  is  most  useful  for  the  detection  of  species 
cither  at  close  range  or  in  high  concentration  in  the  atmosphere  such 
as  N2,  O2,  and  H2O  It  has  a  range  capability  of  tens  of  kilometers 
for  the  detection  of  atmospheric  N2  but  is  limited  to  about  a  few 
hundred  meters  for  pollution  monitoring.  In  addition,  its  sensitivirv 
is  low  (on  the  order  of  100  ppm)  as  compared  with  other 
techniques.  Despite  these  shortcomings,  Raman  scattering  has 
several  features  that  have  made  it  a  useful  remote  sensing  tool.  Its 
most  noteworthy  feature  is  that  the  laser  wavelength  docs  not  have 
to  be  tuned  across  an  absorption  line  since  the  spectral  information 
is  given  by  the  frequency  shiff  of  the  emission,  which  is  independent 
of  the  laser  wavelength.  The  strength  of  the  scattering,  except  in 
certain  resonance  regions,  is  only  weakly  dependent  on  laser  wave¬ 
length,  and  therefore  Raman  lidar  docs  not  require  a  particular  laser 
frequency  as  long  as  the  one  chosen  is  effectively  free  of  atmospheric 
interference. 

The  Raman-scattcrcd  wavclcngth-shiffcd  spectra  arc  distinct  for 
each  molecular  species,  which  allows  separate  analysis  of  the  Raman 
returns  from  each  species.  This  spectral  discrimination  can  be  used 
to  advantage  by  comparing  the  relative  intensity  of  the  returns  from 
the  molecule  that  is  being  dctcacd  to  the  intensity  from  atmospheric 
N2,  whose  concentration  is  known.  This  ratio  approach  simplifies  a 
Raman  lidar  measurement  by  eliminating  atmospheric  or  instru¬ 
mental  uncertainties  in  the  lidar  equation. 

A  wide  variety  of  Raman  lidar  measurements  of  atmospheric 
polluting  species  have  been  made  in  regions  where  their  concentra¬ 
tions  arc  quite  high.  Reviews  of  such  measurements  arc  given  by 
Inaba  and  Kobavashi  (32),  who  discuss  Raman  lidar  measurements 
of  SO2,  NO,  CO,  H2S,  C2H4,  CH4,  and  H2CO  present  in  high 
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concentration  (100  to  1000  ppm)  in  oil  smoke  plumes  and  automo¬ 
bile  exhausts,  at  ranges  on  the  order  of  30  to  100  m.  Another 
application  of  Raman  lidar  arises  from  the  fact  that  the  backscattcrcd 
Raman  intensity  for  a  given  rotational  band  is  strongly  temperature 
dependent.  This  property  has  been  used  for  the  remote  measure¬ 
ment  of  atmospheric  temperature  profiles  at  ranges  of  over  2  km 
with  an  uncertainty  of  less  than  ±2  K.  These  results  are  described  in 
a  review  by  Cooney  of  the  use  of  Raman  scattering  for  the  remote 
sensing  of  meteorological  properties  (33).  In  recent  Raman  lidar 
experiments,  researchers  have  also  studied  water  vapor  profiles  near 
the  water- atmosphere  boundary  layer,  subsurface  ocean  temperature 
profiles,  and  ocean  turbidity.  Long-range  (30  km)  measurements  of 
atmospheric  temperature  by  means  of  Raman  lidar  have  also  been 
reported  (15,  34). 

Doppler  lidar.  Doppler  shifts  in  the  return  lidar  signals  have  been 
used  to  measure  wind  velocities  and  to  differentiate  between  molec¬ 
ular  and  aerosol  returns  in  the  atmosphere.  Although  these  shifts  arc 
small  (a  fractional  change  in  frequency  of  approximately  10“ 8  for  a 
velocity  of  1  m/scc  at  a  wavelength  of  10  pm),  they  can  be  measured 
by  optical  heterodyne  detection  techniques.  Heterodyne  detection 
involves  the  optical  mixing  of  the  lidar  return  with  another  laser 
operating  at  or  near  the  lidar  transmitter  wavelength,  with  detection 
of  the  difference  or  beat  frequency  of  the  mixed  signal. 

Most  Doppler  lidar  measurements  have  used  10- pm  C02  lasers 
because  C02  lasers  offer  high  power  and  stable  single-frequency 
operation  (35-36).  A  review  of  work  in  this  area,  including  a  large 
bibliography,  has  been  given  by  Bilbro  (35)  R:ccndv,  Hall  and  co¬ 
workers  at  NOAA  used  a  C02  Doppler  lidar  to  measure  the  velocity 
and  direction  of  atmospheric  wind  fields  at  ranges  up  to  12  km  (37), 
This  system  has  yielded  information  on  boundary  layer  flow  near 
storm  gust  fronts  and  wind  shears  near  airports;  Fig.  8  shows  a  plot 
of  wind  speed  and  direction  measured  with  Doppler  lidar  near 
Dcnver-Stapleton  Airport  during  a  storm,  showing  the  presence  of  a 
wind  shear-micro  burst  (38).  Doppler  lidar  systems  have  also  been 
used  to  measure  aircraft  vortices  and  clear  air  turbulence.  In 
particular,  an  airborne  C02  laser  system  based  at  the  NASA 
Marshall  center  has  been  used  to  construct  two-dimensional  vector 
wind  fields  in  the  vicinity  of  the  aircraft  in  flight  (39). 
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Fig.  8.  Doppler  lidar  measurement  of  wind  direction  and  velocity  near 
Dcnver-Stapleton  Airport  during  a  storm;  the  arrow  points  to  the  presence 
of  a  strong,  localized  downburst-gustfront.  [Adapted  from  (38) J 
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Doppler  broadening  effects  have  also  been  used  by  a  University  of 
Wisconsin  group  (40)  to  separate  backscattcrcd  lidar  signals  into 
molecular  and  aerosol  components.  This  distinction  can  be  made 
because  the  Doppler  broadening  due  to  the  thermal  motion  of 
molecules  is  about  two  orders  of  magnitude  larger  than  that  of  the 
heavier  aerosols.  The  separation  is  achieved  by  using  a  high- spectral- 
resolution  Fabry-Perot  optical  interferometer  to  separate  the  Dop¬ 
pler-shifted  molecular  and  nearly  unshifted  aerosol  lidar  returns.  The 
system  has  been  used  in  an  aircraft  to  measure  optical  depth  profiles 
and  separate  values  of  the  molecular  and  aerosol  backscatter  cross 
sections  as  a  function  of  height  to  an  altitude  of  1.6  km. 


Accuracy  Limitations 

The  accuracy  of  a  lidar  measurement  is  limited  to  the  accuracy 
with  which  the  average  value  of  the  lidar  returns  can  be  determined, 
and  is  therefore  influenced  by  the  statistical  fluctuations  of  the  lidar 
signals.  These  fluctuations  are  caused  by  several  factors,  including 
detector  noise,  changes  in  the  propagation  and  attenuation  charac¬ 
teristics  due  to  atmospheric  turbulence,  changes  in  the  geometry  of 
the  lidar  beam  and  target,  and  variability  in  the  laser  power  and 
wavelength.  The  limitations  imposed  by  detector  noise  and  signal 
fluctuations  appear  to  be  the  dominant  factors  determining  lidar 
measurement  accuracy  (20). 

The  noise  level  of  the  optical  detector  determines,  in  part,  the 
accuracy  or  sensitivity  of  the  lidar  system;  the  lidar  return  signal 
must  be  greater  than  this  noise  level  to  be  measured  with  precision. 
The  ratio  of  the  lidar  signal  strength  and  the  noise  of  the  system  mav 
be  expressed  analytically  by  noting  that  for  a  direct-detection  lidar 
with  a  photodiode  used  tor  detection  of  the  optical  returns,  the 
detector  voltage  signal- to- noise  (SIN)  ratio  is  given  by  (5,  41) 

S'N  “  [D(Pr  +  PB)  +  fo  +  /e]1'2  (6) 

where  P b  is  the  noise  associated  with  the  background  radiation 
power  incident  on  the  detector  Io  is  the  dark  current  of  the 
detector,  /e  is  the  equivalent  noise  current  of  the  detector  preampli¬ 
fier,  and  C  and  D  arc  constants  related  to  the  signal  bandwidth, 
optical  wavelength,  and  detector  quantum  efficiency.  A  similar 
cquauon  may  be  written  for  heterodyne  detection,  modified  to 
include  a  local  oscillator  power  term,  and  for  optical  detection  with  a 
photomultiplier  tube  or  avalanche  photodiode  to  include  internal 
gain  (6,  41).  The  denominator  in  Eq.  6  is  associated  with  the  shot 
noise  of  the  system  and  is  due  to  quantum  fluctuations  of  the  signal 
and  noise  currents.  As  expected,  the  SIN  ratio  for  different  types  of 
lidar  systems  will  depend  on  the  noise  characteristics  of  the  optical 
detector  and  on  the  spectral  filtering  used  in  the  system.  It  is  often 
possible  to  choose  a  detector  so  that  only  one  noise  term  in  Eq.  6  is 
dominant,  in  this  case,  one  can  classify  the  detection  process  as 
being  limited  by  the  individual  dominant  noise  source,  such  as  signal 
shot  noise,  background  noise,  or  amplifier  noise.  Optical  detection 
in  the  UV,  visible,  and  near  IR  (wavelengths  up  to  1  pm)  uses 
photomultiplier  tubes  that  have,  under  optimal  conditions,  excep¬ 
tionally  low  noise  levels,  are  capable  of  detecting  a  single  photon  of 
fight,  and  arc  essentially  shot  noise-limited.  Optical  detection  in  the 
middle  IR  (1  to  10  pm)  usually  uses  semiconductor  materials  for 
direct  detection  of  the  optical  radiation;  these  devices  are  often 
background  noise-limited  and  have  noise  levels  that  arc  several 
orders  of  magnitude  higher  than  photomultiplier  tube  noise  levels. 
Enhanced  detection  in  the  middle  IR  can  be  accomplished  with 
optical  heterodyne  detection,  which  is  capable  of  single  photon 
detection,  but  only  with  a  considerable  increase  in  lidar  complexity 
due  to  the  addition  of  another  laser  used  for  optical  mixing  with  the 
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lidar  returns  and,  in  the  ease  of  a  DIAL  system,  the  precision 
tunabilitv  and  tracking  of  these  laser  wavelengths. 

The  accuracy  of  a  lidar  measurement  is  related  to  the  lidar  S/N 
ratio  shown  in  Eq.  6  as  ct0  =  (S/N)"1,  where  cr0  is  the  normalized 
statistical  standard  deviation  of  the  fluctuating  lidar  returns  and  is  a 
measure  of  the  uncertainty  in  the  estimate  of  the  average  value  of  the 
returns  (41).  It  is  well  known  that  the  accuracy  of  a  measurement 
can  be  enhanced  by  signal  averaging  and  that  the  S/N  ratio  or 
accuracy  will  increase  as  the  square  root  of  the  number  of  pulses 
integrated.  This  may  be  expressed  analytically  as  an  =  oV»1/2,  where 
c 7„  is  the  standard  deviation  in  the  estimate  of  the  mean  value  of  the 
n  pulses.  This  w17  improvement  in  accuracy  by  signal  averaging  can 
be  anticipated  only  if  the  measured  returns  from  successive  signals 
are  uncorrelatcd,  which  means  that  the  statistical  properties  of  all  the 
terms  in  Eq,  6  be  temporally  random  (stochastic).  In  many  lidar 
examples,  the  dominant  noise  sources  are  random  and  an  «l/2 
improvement  is  observed.  Such  lidar  measurements  include,  as 
examples,  Raman  and  fluorescence  lidar  measurements  of  the  upper 
atmosphere,  DIAL  measurements  when  atmospheric  attenuation 
fluctuations  are  negligible,  and  the  reduction  of  speckle  (interfer¬ 
ence)  effects  in  Doppler  lidar  measurements. 

In  some  lidar  measurements,  however,  the  fluctuations  observed 
are  not  random  and  arc  caused  bv  long-term  variability  in  the  optical 
properties  and  attenuation  of  the  atmosphere  during  the  measure¬ 
ment  period.  This  is  often  seen  in  DIAL  measurements  of  the 
atmosphere,  especially  when  the  laser  beam  is  transmitted  in  the 
lower  atmosphere  so  that  it  is  most  susceptible  to  the  influence  of 
long-term  drifts  in  the  composition  and  temperature  of  the  atmo¬ 
sphere  near  the  surface  of  the  earth.  Under  these  conditions  the 
statistical  distribution  and  temporal  characteristics  of  the  terms  in 
Eq.  6  arc  subject  to  both  random  and  nonrandom  influences,  and  a 
deviation  from  the  expected  nV2  improvement  in  the  S/N  ratio  or  ct„ 
is  observed  (42).  An  example  of  this  effect  is  presented  in  Fig.  9, 
which  shows  a  logarithmic  plot  of  ct„  as  a  function  of  «,  the  number 
of  lidar  pulses  averaged,  for  DIAL  data  obtained  with  a  CO2  laser; 
for  comparison,  a  rv2  improvement  is  shown  as  the  straight  line  in 
the  figure.  The  marked  deviation  from  nl/ 2  behavior  shown  in  Fig,  9 
indicates  the  extent  to  which  nonrandom  processes  may  reduce  the 
effectiveness  of  signal  averaging.  Figure  9  is  only  one  example  of  this 
observed  phenomenon,  in  fact,  the  deviation  from  nm  behavior 
ranges  from  zero  to  values  significantly  greater  than  that  shown  The 
physical  sources  of  the  correlation  are  believed  to  be  atmospheric 
turbulence  and  long-term  drifts,  which  lead  to  temporal  changes  in 
atmosphere  attenuation  during  the  measurement  interval.  Both  of 
these  effects  would  produce  nonrandom  fluctuations  in  the  atmo¬ 
spheric  attenuation,  a,  and  thus  influence  the  statistical  properties  of 
Pr 

In  order  to  better  understand  the  data  shown  in  Fig.  9,  we  made  a 
theoretical  analysis  to  quantify  the  effect  of  correlation  on  a  statisti¬ 
cal  ensemble  of  data,  with  direct  application  to  DIAL  (43).  This 
analysis  showed  that  the  inclusion  of  correlation  changes  the  simple 
w"1/2  dependence  of  <j„  to 

=  3?!1  +  2  2(1  (7) 

n  j-i 

where  p;  is  the  temporal  autocorrelation  coefficient  for  a  delay  time 
Jt,  and  t  is  the  time  interval  between  successive  lidar  measurements. 
For  stochastic  data,  pj  is  zero  and  Eq.  7  reduces  to  the  n~ 1/2  form  as 
expected.  As  seen  in  Eq.  7,  the  effect  of  a  nonzero  value  of  pj  is  to 
reduce  the  improvement  in  signal  averaging.  Further  DIAL  experi¬ 
ments  (43)  have  shown  that  Eq.  7  is  in  excellent  agreement  with  the 
observed  DIAL  data  as  presented  in  Fig.  9. 

Equation  7  can  be  shown  to  reduce  to  cr„  —  a0  p>1/2  for  large  n  and 
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Fig.  9.  Measurement  ac¬ 
curacy  (percent  standard 
deviation)  of  C02  lidar 
signal  returns  from  a  dif¬ 
fusely  reflecting  target  at 
a  range  of  2.7  km  as  a 
function  of  the  number 
of  pulses  averaged.  The 
reduced  accuracy  im¬ 
provement  relative  to 
n~xa  behavior  is  due  to 
long-term  drifts  in  atmo¬ 
spheric  attenuation. 

nonzero  pj.  This  result  shows  that  a  limit  exists,  independent  of  w, 
for  the  improvement  in  the  accuracy  of  a  lidar  measurement  if 
nonrandom  changes  occur  in  the  atmosphere  during  the  measure¬ 
ment  interval  (43).  Equation  7  is  valid,  in  general,  for  all  experimen¬ 
tal  measurements  and  describes  quantitatively  the  often-observed 
diminished  improvement  from  long-term  averaging  when  nonran¬ 
dom  perturbations  are  causing  changes  in  the  average  signal  during 
the  measurement  period. 

Finally,  while  Eq  7  docs  predict  the  increase  in  accuracy  of  a  lidar 
measurement  through  signal  averaging,  <7„  can  be  reduced  by 
techniques  that  reduce  either  the  initial  values  of  <j0  or  the  initial 
value  of  the  correlation  coefficient.  These  techniques  (43—45) 
include  the  use  of  very  high  repetition  rate  lasers  to  increase  the 
number  of  lidar  samples  collected  before  the  atmosphere  changes 
appreciably;  the  use  of  dual,  nearly  simultaneously  transmitting, 
lasers  for  DIAL  experiments,  and  the  use  of  laser  wavelength  agility. 
Even  with  some  of  these  techniques,  the  accuracy  of  most  current 
lidar  experiments  has  been  limited  to  values  on  the  order  of  a  percent 
(that  is,  (7„  *»  1%)  under  optimal  conditions  and  to  values  on  the 
order  of  1  to  10%  under  normal  atmospheric  conditions 


Future  Directions 

Laser  remote  sensing  is  a  useful  tool  for  probing  the  composition 
and  physical  state  of  the  atmosphere  and,  by  means  of  single-ended 
measurements  at  a  distance,  permits  measurements  at  locations  that 
cannot  be  reached  with  conventional  detection  methods  such  as 
chemical  or  optical  point  sensors  Although  technical  limitations  do 
exist,  in  many  instances,  lidar  offers  the  only  realistic  measurement 
technique  that  can  detect  a  physical  quantity  and  produce  a  three- 
dimensional  map  of  the  quantity  in  real  time. 

The  future  of  laser  remote  sensing  appears  to  be  quite  promising 
and  is  dependent  upon  factors  that  include  (i)  the  development  of 
practical,  new,  eye-safe  laser  sources  that  cover  the  spectral  gaps 
shown  in  Fig.  4,  (ii)  further  reduction  of  lidar  system  complexity, 
cost,  and  size,  and  (iii)  the  future  use  of  lidar  techniques  in  new 
applications  Among  the  latter  arc  several  emerging  industrial  and 
medical  applications  that  include  the  use  of  fiber  optics  in  enclosed 
or  hostile  environments,  the  use  of  airborne  lidars  for  oil,  gas,  and 
geological  exploration,  and  the  exploitation  of  the  long-range 
(>100  km)  detection  capability  and  wide  surveillance  potential  of 
space-borne  laser  remote  sensing.  Examples  that  are  currently  being 
explored  include  (i)  the  detection  of  methane  gas  leaks  in  a  coal  mine 
with  a  diode  laser  lidar  system,  (ii)  the  use  of  a  laser  coupled  to  a 
low-loss  optical  fiber  network  to  detea  methane  or  natural  gas  leaks 
in  an  industrial  processing  plant,  (iii)  the  Doppler  lidar  measure¬ 
ment  of  global  wind  fields  from  a  satellite  to  increase  the  accuracy  of 
weather  forecasts,  and  (iv)  the  planned  use  of  lidar  aboard  the 
NASA  space-borne  Earth  Observing  System  for  global  temperature, 
water  vapor,  and  pressure  measurements. 
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